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Gravitational microlensing occurs when a foreground star happens to pass very close to our line 
of sight to a more distant background star. The foreground star acts as a lens, splitting the light from 
the background source star into two images, which are typically unresolved. However, these images 
of the source are also magnified, by an amount that depends on the angular separation between the 
lens and source. The relative motion between the lens and source therefore results in a time- variable 
magnification of the source: a microlensing event. If the foreground star happens to host a planet with 
projected separation near the paths of these images, the planet will also act as a lens, further perturbing 
the images and resulting in a characteristic, short-lived signature of the planet. This chapter provides 
an introduction to the discovery and characterization of exoplanets with gravitational microlensing. 
The theoretical foundation of the method is reviewed, focusing in particular on the phenomenology of 
planetary microlensing perturbations. The strengths and weaknesses of the microlensing technique are 
discussed, highlighting the fact that it is sensitive to low-mass planetary companions to stars throughout 
the Galactic disk and foreground bulge, and that its sensitivity peaks for planet separations just beyond the 
snow line. An overview of the practice of microlensing planet searches is given, with a discussion of some 
of the challenges with detecting and analyzing planetary perturbations. The chapter concludes with a 
review of the results that have been obtained to date, and a discussion of the near and long-term prospects 
for microlensing planet surveys. Ultimately, microlensing is potentially sensitive to multiple-planet 
systems containing analogs of all the solar system planets except Mercury, as well as to free floating 
planets, and will provide a crucial test of planet formation theories by determining the demographics of 
planets throughout the Galaxy. 



1. INTRODUCTION 

Gravitational lensing generally refers to the bending of 
light rays of a background light source by a foreground 
mass. Gravitational microlensing, on the other hand, tra- 
ditionally refers to the special case when multiple images 
are created but have separations of less than a few milliarc- 
seconds, and hence are unresolved with current capabilities. 
Although the idea of the gravitational deflection of light by 
massive bodies well predates General Relativity, and can 
be traced as far back as Sir Isaac NewtorQ, the concept of 
gravitational microlensing appears to be attributable to Ein- 
stein himself. In 1936 Einstein published a paper in which 
he derived the equations of microlensing by a foreground 
star closely aligned to a background star {Einstein, 1936). 
Indeed it seems that Einstein had been thinking about this 
idea as far back as 1912 (Renn et ai, 1997), and perhaps 
had even hoped to use the phenomenon to explain the ap- 
pearance of Nova Geminorum 1912 (Sauer, 2008). How- 
ever, by 1936 he had dismissed the practical significance of 
the microlensing effect, concluding that "there is no great 
chance of observing this phenomenon" (Einstein, 1936). 

Indeed, there is no "great chance" of observing gravita- 
tional microlensing. The optical depth to gravitational mi- 
crolensing, i.e., the probability that any given star is being 
appreciably lensed at any given time, is of order 10^^ to- 



See Schneider et al. 1992 for a thorough recounting of the history of grav- 
itational lensing. 



ward the Galactic bulge, and is generally similar or smaller 
for other lines of sighjl Thus at least partly due to this 
low probability, the idea of gravitational microlensing lay 
mostly dormant for five decades after Einstein's 1936 pa- 
per (with some notable exceptions, e.g., Liebes, 1964; Refs- 
dal, 1964). It was not until the seminal paper by Paczynski 
(1986) that the idea of gravitational microlensing was resur- 
rected, and then finally put into practice with the initiation 
of several observational searches for microlensing events 
toward the Large and Small Magellenic Clouds and Galac- 
tic bulge in the early 1990s (Alcock et al. 1993; Aubourg 
et al. 1993; Udalski et al, 1993). The roster of detected 
microlensing events now numbers in the thousands. These 
events have been discovered by many different collabora- 
tions, toward several lines of sight including the Magellenic 
Clouds (Alcock et al., 1991, 2000; Palanque-Delabrouille 
et al., 1998), and M31 (Paulin-Henriksson et al. 2002; de 
Jong et al., 2004; Uglesich et al., 2004; Calchi Novati et 
al., 2005), with the vast majority found toward the Galactic 
bulge (Udalski et al, 2000; Sumi et al., 2003; Thomas et 
al., 2005; Hamadache et al., 2006) or nearby fields in the 
Galactic plane (Derue et al., 2000). 



^The phenomenon of gravitational lensing of multiply-imaged quasars by 
stars in the foreground galaxy that is creating the multiple images of the 
quasar is also referred to as microlensing. In this instance the optical depth 
to microlensing can be of order unity. However, in this chapter we are 
concerned only with gravitational microlensing of star's within our Galaxy 
or in nearby galaxies, where the optical depths to microlensing are always 
small. 
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Fig. 1. — Left: The images (dotted ovals) are shown for several different positions of the source (solid circles), along with the primary 
lens (dot) and Einstein ring (long dashed circle). If the primary lens has a planet near the path of one of the images, i.e. within the 
short-dashed lines, then the planet will perturb the light from the source, creating a deviation to the single lens light curve. Right: The 
magnification as a function of time is shown for the case of a single lens (solid) and accompanying planet (dotted) located at the position 
of the X in the left panel. If the planet was located at the + instead, then there would be no detectable perturbation, and the resulting 
light curve would be essentially identical to the solid curve. 



Even before the first microlensing events were detected, 
it was suggested by Mao & Paczynski (1991) that gravita- 
tional microlensing might also be used to discover plane- 
tary companions to the primary microlens stars. The basic 
idea is illustrated in Figure[T] As the foreground star passes 
close to the line of sight to the background source, it splits 
the source into two images, which sweep out curved trajec- 
tories on the sky as the foreground lens star passes close to 
the hne of sight to the source. These images typically have 
separations of order one milliarcsecond, and so are unre- 
solved. However, the total area of these images is also larger 
than the area of the source, and as a result the background 
star also exhibits a time-variable magnification, which is re- 
feiTed to as a microlensing event. If the foreground star 
happens to host a planet with projected separation near the 
paths of one of these two images, the planet will further 
perturb the image, resulting in a characteristic, short-lived 
signature of the planet. 

Gould & Loeb (1992) considered this novel method of 
detecting planets in detail, and laid out the practical re- 
quirements for an observational search for planets with 
microlensing. In particular, they advocated a "two tier" 
approach. First, microlensing events are discovered and 
alerted real-time by a single survey telescope that monitors 
many square degrees of the Galactic bulge on a roughly 
nightly basis. Second, these ongoing events are then 
densely monitored with many smaller telescopes to dis- 



cover the short-lived signatures of planetary perturbations. 

The search for planets with microlensing began in 
earnest in 1995 with the formation of several follow-up 
collaborations dedicated to searching for planetary devia- 
tions in ongoing events {Alb row et ai, 1998; Rhie et al. 
2000). These initial searches adopted the basic approach 
advocated by Gould & Loeb (1993): monitoring ongo- 
ing events alerted by several survey collaborations (e.g., 
Udalski et al., 1994; Alcock et al. 1996), using networks 
of small telescopes spread throughout the southern hemi- 
sphere. Although microlensing planet searches have ma- 
tured considerably since their initiation, this basic approach 
is still used to this day, with the important modification 
that current follow-up collaborations tend to focus on high- 
magnification events, which individually have higher sensi- 
tivity to planets {Griest & Safizadeh, 1998), as discussed in 
detail below. 

From 1995-2001, no convincing planet detections were 
made, primarily because the relatively small number of 
events being alerted each year (^ 50 — 100) by the sur- 
vey collaborations meant that there were only a few events 
ongoing at any given time, and often these were poorly 
suited for follow-up. Although interesting upper limits were 
placed on the frequency of Jovian planets (Gaudi et al. 
2002, Snodgrass et al., 2004), perhaps the most important 
result during this period was the development of both the 
theory and practice of the microlensing method, which re- 
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suited in its transformation from a theoretical abstraction to 
a viable, practical method of searching for planets. 

In 2001, the OGLE collaboration (Udalski, 2003) up- 
graded to a new camera with a 16 times larger field of view 
and so were able to monitor a larger area of the bulge with 
a higher cadence. As a result, in 2002 OGLE began alert- 
ing nearly an order of magnitude more events per year than 
previous to the upgrade. These improvements in the alert 
rate and cadence, combined with improved cooperation and 
coordination between the survey and follow-up collabora- 
tions, led to the first discovery of an extrasolar planet with 
microlensing in 2003 (Bond et al. 2004). The MOA col- 
laboration upgraded to a 1.8m telescope and 2 deg'^ cam- 
era in 2004 (Sako et al, 2008), and in 2007 the OGLE and 
MOA collaborations started alerting ^ 850 events per year, 
thus ushering in the "golden age" of microlensing planet 
searches. 

The first detections of exoplanets with microlensing 
(Bond et al., 2004; Udalski et al., 2005; Beaulieu et al., 
2006; Gould et al., 2006; Gaudi et al., 2008a; Bennett et 
al., 2008; Dong et al., 2009; Janczak et al., 2009; Sumi 
et al, 2009) have provided important lessons about the 
kinds of information that can be extracted from observed 
events. When microlensing planet surveys were first be- 
ing developed, it was thought that the primary virtue of the 
method would be solely its ability to provide statistics on 
large-separation and low-mass planets. Individual detec- 
tions would be of little interest because the only routinely 
measured property of the planets would be the planet/star 
mass ratio, and the nature of the host star in any given 
system would remain unknown because the microlensing 
event itself provides little information about the properties 
of the host star, and follow-up reconnaissance would be 
difficult or impossible due to the large distances of the de- 
tected systems from Earth. In fact, experience with actual 
events has shown that much more information can typically 
be gleaned from a combination of a detailed analysis of the 
light curve, and follow-up, high-resolution imaging. As a 
result, in most cases it is possible to infer the mass of the 
primary lens and the basic physical properties of planetary 
system, including the planet mass and physical separation. 
In some special cases it is possible to infer considerably 
more information. 

The theoretical foundation of microlensing is highly 
technical, the practical challenges associated with search- 
ing for planets with microlensing are great, and the analysis 
of observed microlensing datasets is complicated and time- 
consuming. Furthermore, although it is possible to learn far 
more about individual systems than was originally thought, 
extracting this information is difficult, and it is certainly the 
case that the systems will never be as well studied as planets 
systems found by other methods in the solar neighborhood. 

Give the intrinsic limitations and difficulties of the mi- 
crolensing method, and might therefore wonder: "why 
bother?" The primary motivation for going to all this trou- 
ble is that microlensing is sensitive to planets in a region of 
parameter space that is difficult or impossible to probe with 



other methods, namely low-mass planets beyond the 'snow 
line' , the point in the protoplanetary disk beyond which ices 
can exist. The snow line plays a crucial role in the currently- 
favored model of planet formation. Thus, even the first 
few microlensing planet detections have provided impor- 
tant constraints on planet formation theories. The second 
generation of microlensing surveys will potentially be sen- 
sitive to multiple-planet systems containing analogs of all 
the solar system planets except Mercury, as well as to free 
floating planets. Ultimately, when combined with the re- 
sults from other complementary surveys, microlensing sur- 
veys can potentially yield a complete picture of the demo- 
graphics of essentially all planets with masses greater than 
that of Mars. Thus it is well worth the effort, even given the 
drawbacks. 

The primary goal of this chapter is to provide a general 
introduction to gravitational microlensing and searches for 
planets using this method. Currently, the single biggest ob- 
stacle to the progress of microlensing searches for planets 
is simply a lack of human power There are already more 
observed binary and planetary events than can be modeled 
by the handful of researchers in the world with sufficient 
expertise to do so. This situation is likely to get worse 
when next-generation microlensing planet surveys come on 
line. However, because the theoretical foundation of mi- 
crolensing is quite technical and the practical implementa- 
tion of microlensing planet searches is fairly complex, it can 
be very difficult and time consuming for the uninitiated to 
gain sufficient familiarity with the field to make meaning- 
ful contributions. This difficulty is exacerbated by the fact 
that the explanations of the relevant concepts are scattered 
over dozens of journal articles. The aim of this chapter is 
to partially remedy this situation by providing a reasonably 
comprehensive review, which a beginner can use to gain a 
least a basic familiarity with the many concepts, tools, and 
techniques that are required to actively participate and con- 
tribute to current microlensing planet searches. Readers are 
encouraged to peruse reviews by Paczynski 1996, Sackett 
1999, Mao 2008, and Bennett 2009a for additional back- 
ground material. 

SectionlUprovides an overview of the theory of gravita- 
tional microlensing searches for planets. This is the heart 
of the chapter, and is fairly long and detailed. Much of this 
section may not be of interest to all readers, and some of 
the material may be too technical for a beginner to the sub- 
ject. Therefore, the first subsection ( 12.11 ) provides a primer 
on the basic properties and features of microlensing. Be- 
ginners, or casual readers who are only interested in a basic 
introduction to the method itself and its features, can read 
only this section (paying particular attention to Figures[Tl|2l 
andO, and then skip to Section |4] Section [3] discusses the 
practical implementation of planetary microlensing, while 
Section |4] reviews the basic advantages and drawbacks of 
the method. Section |5] provides a summary of the results 
to date, as well as a brief discussion of their implications. 
Section|6]discusses future prospects for microlensing planet 
searches, and in particular the expected yields of a next- 
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Fig. 2. — Left: The lens (L) at a distance D; from tiie observer (O) deflects ligiit from the source (S) at distance Da by the Einstein 
bending angle aa.- The angular positions of the images Q and unlensed source /3 are related by the lens equation, /3 = Q oca, ~ 
— {Ds — Di) / Ds&d- For a point lens, ad ~ 4,GM / {(? DiO). Right: Relation of higher-order observables, the angular (Se) and 
projected {f^) Einstein radii, to physical characteristics of the lensing system. Adapted from Gould, 2000. 



generation ground-based planet search, and a space-based 



2. FOUNDATIONAL CONCEPTS AND EQUATIONS 
2.1. Basic Microlensing 

This section provides a general overview of the basic 
equations, scales, and phenomenology of microlensing by 
a point mass, and a brief introduction to how microlensing 
can be used to find planets, and how such planet searches 
work in practice. It is meant to be self-contained, and there- 
fore the casual reader who is not interested in a detailed dis- 
course on the theory, phenomenology, and practice of plan- 
etary microlensing can simply read this section and then 
skip to Section|4]without significant loss of continuity. 

A microlensing event occurs when a foreground "lens" 
happen to pass very close to our line of sight to a more 
distant background "source." Microlensing is a relatively 
improbable phenomenon, and so in order to maximize the 
event rate, microlensing survey are typically carried about 
toward dense stellar fields. In particular, the majority of mi- 
crolensing planet surveys are carried out toward the Galac- 
tic center. Therefore, for our purposes, the lens is typi- 
cally a main-sequence star or stellar remnant in the fore- 
ground Galactic disk or bulge, whereas the source is a main- 
sequence star or giant typically in the bulge. 

The left panel of Figure |2] shows the basic geometry of 
microlensing. Light from the source at a distance Ds is 



deflected by an angle &d by the lens at a distance Di. For 
a point lens, ad = AGM /{c^DiO), where M is the mass 
of the lens, and 6 is the angular separation of the images of 
the source and the lens on the sk}|j. The relation between 
0, and the angular separation /? between the lens and source 
in the absence of lensing, is called the lens equation, and 
is given trivially by /3 = 6* — ad- From basic geometry 
and using the small-angle approximation, ad{Ds — Di) = 
adDs- Therefore, for a point lens. 



4GM D,-D, 



^6 DsD, 



(1) 



The left panel of Figure[3]shows the basic source and im- 
age configurations for microlensing by a single point mass. 
From Equation [U if the lens is exactly aligned with the 
source (/3 = 0), it images the source into an "Einstein ring" 
with a radius 9^ = ^/kMtt^, where M is the mass of the 
lens, TTici = AU/Drei is the relative lens-source parallax, 
Z^ici = {Dj^^ — Dj^)~^ is the relative lens-source distance, 
and n = 4G/c^AU ~ 8.14 mas Mq^. It is also instructive 



•^This form for the bending angle can be derived heuristically by assuming 
that a photon passing by an object of mass M at a distance b = DiO 
will experience an impulse given by the Newtonian acceleration GM/b^ 
over a time 2fe/c, thereby inducing a velocity perpendicular to the original 
teajectory of 5v = {GM /b'^){2b/ c) = 2GM/{bc). The deflection is 
then &v/c = 2GM/ (bc^). The additional factor of two cannot be derived 
classically, and arises from General Relativity (see, e.g., Schneider et al. 
1992). 
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Fig. 3. — Basic point-mass microlensing. (Left) All angles are normalized by the angular Einstein ring radius 6e, shown as a dashed 
circle of radius 6e- The source (S) is located at an angular separation of u = 0.2 from the lens (L). Two images are created, one image 
outside the Einstein ring (I+), on the same side of the lens as the source with position from the lens of y+ = 0.5(\/ + 4 + it), and 
one inside the Einstein ring, on the opposite side of the lens as the source with position from the lens of y_ = — 0.5(Vw^ + 4 — u). The 
images are compressed radially but elongated tangentially. Since surface brightness is conserved, the magnification of each image is just 
the ratio of its area to the area of the source. Since the images are typically unresolved, only the total magnification of the two images is 
measured, which depends only on u. (Right) Magnification as a function of time (light curves), for the ten trajectories shown in the left 
panel with impact parameters uo ~ 0.01, 0.1, 0.2, 1.0. Time is relative to the time to of the peak of the event (when u = wo), and in 
units of the angular Einstein crossing time ts- Higher magnification implies more elongated images, which leads to increased sensitivity 
to planetary companions. Adapted from Paczynski 1996. 



to note that = ^/2R^~ijD^i, where Rsch = ^GM jc? 
is the Schwarzschild radius of the lens. Quantitatively, 
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which corresponds to a physical Einstein ring radius at the 
distance of the lens of. 
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ways inside the Einstein ringj. As can be seen in Figure [3] 
the angular separation between these images at the time of 
closest alignment is is ^ 26'e- Thus for typical lens masses 
(0.1 — 1 Mq) and lens and source distances (1 — 10 kpc), 
^E ^ mas and so the images are not resolved. However, 
the images are also distorted, and since surface brightness 
is conserved, this implies they are also magnified. The mag- 
nification of each image is just the ratio of the area of the 
image to the area of the source. As can be seen in Fig- 
iJ ure[3] the images are typically elongated tangentially by an 
amount y±/u, but are compressed radially by an amount 
dz/i/du. The magnifications are then 
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(4) 



where x = Di/Dg- 

Normalizing all angles on the sky by 6'e, we can define 
u = /3/6'e and y — 9/6^- Using these definitions and the 
definition for 6'e, the lens equation ^ reduces to the form 

u^y-y~^- (5) 

This is equivalent to a quadratic equation in y: y'^—uy— 1 — 
0. Thus in the case of imperfect alignment (u ^ 0), there 
are two images, with positions 
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(6) 



The positive, or "major" image is always outside the Ein- 
stein ring, whereas the negative, or "minor" image is al- 



*The images formed by a gravitational lens can also be found by determin- 
ing the relative time delay as a function of the (vector) angle for hypo- 
thetical light rays propagating from the source. This time delay function 
includes the effect of the difference in the geometric path length, as well 
the gravitational (Shapiro) delay, as a function of 0. The images are then 
located at the stationary points (maxima, minima, or saddle points) of the 
time delay surface. For a point-mass lens, the positive solution in Equa- 
tion|6]con'esponds to a minimum in the time delay surface, and so is also 
referred to as the "minimum" image. The negative solution corresponds to 
a saddle point, and so is also refeiTed to as the "saddle" image. There is 
formally a third image, coiresponding to the maximum of the time delay 
surface, but this image is located behind the lens, and is infinitely demag- 
nified (for a true point lens). 
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and thus the total magnification is 



L^/v? + 4 



(8) 



Note that, as u — > oo, — > 1 and yl_ — > 0. Also, for 
7i ^ 1, the magnification takes the form A ~ 1/u. 

Because the lens, source, and observer are all in relative 
motion, the angular separation between the lens and source 
is a function of time. Therefore, the magnification of the 
source is also a function of time: a microlensing event. In 
the simplest case of uniform, rectilinear motion, the relative 
lens-source motion can be parametrized by 



u{t) 



t-tp 



1/2 



(9) 



where uq is the minimum separation between the lens and 
source (in units of 6'e), is the time of maximum mag- 
nification (corresponding to when u = uq), and <e is 
the time scale to cross the angular Einstein ring radius, 
= ^?E//"icb where /i^ci is the proper motion of the source 
relative to the lens. This form for u gives rise to a smooth, 
symmetric microlensing event with a characteristic form 
(often called a "Paczynski curve"), as shown in the right 
panel of Figure |3] The typical timescales for events toward 
the Galactic bulge are of order a month. 
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but can range from less than a day to years. The magnifica- 
tion is A > 1.34 for u < 1, and so the magnifications are 
substantial and easily detectable. 

If the lens star happens to have a planetary companion, 
and the planet happens to be located near the paths of one 
or both of the two images created by the primary lens, the 
companion can create a short-timescale perturbation on the 
primary microlensing event when the image(s) sweep by 
the planet (Mao & Paczynski 1991; Gould & Loeb, 1992). 
There are two conceptually different channels by which 
planets can be detected with microlensing, delineated by 
the maximum magnification Amax of the primary event in 
which the planet is detected. 

Consider the case of a relatively low-magnification (i.e., 
^max ^ lOorwQ > 0.1) primary event, such as illustrated 
in Figure [T] Over the course of the primary event, the two 
images sweep out a path on the sky, with the path of the 
major (minimum) image entirely contained outside the Ein- 
stein ring radius, and the path of the minor (saddle) image 
entirely contained within the Einstein ring radius. In the 
main detection channel, the primary hosts a planet which 
happens to be located near the path of one|f| of these two 
images for such a low-magnification event. As the image 



sweeps by the position of the planet, the planet will further 
perturb the light from this image and yield a short-timescale 
deviation (Gould & Loeb, 1992). The duration of this devia- 
tion is ^ iE,p — q^^^t-E, where q — mp/M is the mass ratio 
and nip is the planet mass, and the magnitude of the pertur- 
bation depends on how close the perturbed image passes 
by the planet. Given a range of primary event durations of 
^ 10—100 days, the duration of the perturbations range 
from a few hours for an Earth-mass planet to a few days 
for Jovian-mass planets. As can be seen in Figure |2] the 
location of the perturbation relative to the peak of the pri- 
mary event depends on two parameters: a, the angle of the 
projected star-planet axis with respect to the source trajec- 
torjO, and d, the instantaneous angular separation between 
the planet and host star in units of Be- 

Because the orientation of the source trajectory relative 
to the planet position is random, the time of this perturba- 
tion is not predictable and the detection probability is ^ 
^(^o,p)^'e.p/^'e, where A{to^p) is the unperturbed magnifi- 
cation of the image that is being perturbed at the time io,p of 
the perturbation (Home et al, 2009), and 6'e,p = q^^^Oe is 
the Einstein ring radius of the planet. Here the factor A ac- 
counts for the fact that the area of the image plane covered 
by magnified images is larger by their magnification. 

Since the planet must be located near one of the two pri- 
mary images in order to yield a detectable deviation, and 
these images are always located near the Einstein ring ra- 
1 dius when the source is significantly magnified (see Fig- 
,ure[T]l, the sensitivity of the microlensing method peaks for 
planet-star projected separations of ^ te, i.e., for d ^ 1. 
However, microlensing can also detect planets well outside 
the Einstein ring (d ^ 1), albeit with less sensitivity. Since 
the magnification of the minor image decreases with posi- 
tion as (see Equations |6] and |7]), microlensing is gener- 
ally not sensitive to planets d 1, i.e., close-in planets. 

Given the existence of a planet with a projected separa- 
tion within a factor of ^ 2 of the Einstein ring radius, the 
detection probabilities range from tens of percent for Jo- 
vian planets to a few percent for Earth-mass planets (Gould 
& Loeb, 1992; Bolatto & Falco, 1994; Bennett & Rhie, 
1996; Peale, 2001). Detecting planets via the main chan- 
nel requires substantial commitment of resources because 
the unpredictable nature of the perturbation requires dense, 
continuous sampling, and furthermore the detection prob- 
ability per event is relatively low so many events must be 
monitored. 

A useful feature of low-magnification planetary mi- 
crolensing events such as that shown in Figure [T] is that it is 
possible to essentially 'read off the light curve parameters 
from the observed features (Gould & Loeb, 1992; Gaudi 
& Gould, 1997b). For the primary event, the three gross 
observable parameters are the time of maximum magnifica- 



^As can be seen in Figure[T] for low-magnification events, tlie two images 
are well-separated, and thus a low mass ratio planet wiU generally only 
significantly perturb one of the two images. 



® Conventionally, both the lensing deflection angle and the source trajectory 
angle are denoted by the symbol a. In order to maintain contact with 
the literature, this unfortunate convention is maintained here, but to avoid 
confusion the deflection angle is denoted with the subscript "d". 
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tion to, the peak magnification Amax, and a measure of the 
duration of the event such as its full-width half-maximum 
tpwHM- The latter two observables can then be related to 
the Einstein timescale <e and impact parameter ua using 
Equations [8] and |9]0. For example, for small uq we have 
uo ^ ^max and O.btpwHMUQ^. Unfortunately, 

only t-E contains any information about the physical proper- 
ties of the lens, and then only in a degenerate combination 
of the lens mass, distance, and relative lens-source proper 
motion. However, as discussed in detail in Section 13.31 in 
many cases it is often possible to obtain additional infor- 
mation which partially or totally breaks this degeneracy. In 
particular, in those cases where it is possible to isolate the 
light from the lens itself, this measurement can be used to 
constrain the lens mass {Bennett et ai, 2007). The three 
parameters that characterize the planetary perturbation are 
the duration of the perturbation, the time of the perturba- 
tion to^p, and the magnitude of the perturbation, Sp. As 
mentioned previously, the duration of the perturbation is 
proportional to t^.p = q^^'^tE, and so gives the planet/star 
mass ratio q. The time and magnitude of the perturbation 
then specify the projected separation d and angle of the 
source trajectory with respect to the binary axis a, since 
to p and 5p depend on the location of the planet relative to 
the path of the perturbed image (Gaudi & Gould, 1997b). 

The other channel by which microlensing can detect 
planets is in high magnification events (Griest & Safizadeh, 
199811 In addition to perturbing images that happen to pass 
nearby, planets will also distort the perfect circular symme- 
try of the Einstein ring. Near the peak of high-magnification 
events, as the lens passes very close to the observer-source 
Une of sight (i.e. when u <C 1), the two primary images 
are highly elongated and sweep along the Einstein ring, 
thus probing this distortion. For very high-magnification 
events, these images probe nearly the entire Einstein ring 
radius and so are sensitive to all planets with separations 
near te, regardless of their orientation with respect to the 
source trajectory. Thus high-magnification events can have 
nearly 100% sensitivity to planets near the Einstein ring 
radius, and are very sensitive to low-mass planets (Griest 
& Safizadeh, 1998). However, these events are rare: a 
fraction ^ l/^max of events have maximum magnifica- 
tion > Amax- Fortunately, these events can often be pre- 
dicted several hours to several days ahead of peak, and fur- 
thermore the times of high sensitivity to planets are within 
a full-width half-maximum of the event peak, or roughly 
a day for typical high-magnification events (Rattenbury et 
at, 2002). Thus scarce observing resources can be con- 
centrated on these few events and only during the times of 
maximum sensitivity. Because the source stars are highly 
magnified, it is also possible to use more common, smaller- 

^For the purposes of exposition, this discussion ignores blending, which is 
generally important and complicates the interpretation of observed light 
curves. See Section im 

*There is no formal definition for 'high'-magnilication events, however typ- 
ically high-magnification refers to events with maximum magnification 
Amax > 100, corresponding to impact parameters uq < 0.01. 



aperture telescopes. 

2.2. Theory of Microlensing 

Gravitational lensing can be thought of as the mapping 
(3^6 between the angular position of a source /3 in the ab- 
sence of lensing to the angular position(s) of the image(s) 
of the source under the action of the gravitational lens. This 
mapping is given by the lens equation, 

f3 = e-adie), (11) 

where a.^ is deflection of the source due to the lens. 

Consider a source at a distance Dg and a lens located at 
a distance Di from the observer Figure|2]shows the lensing 
geometry. The deflection angle cy.d is related to the angle 
&d by which the lens mass bends the light ray from the 
source by &d{Ds — Di) — aaDs- Assume the lens is a 
system of Nl point masses each with mass mi and position 
9,n.i- Further assume that the lenses are static (or, more 
precisely, moving much more slowly than c), and their dis- 
tribution along the line of sight is small in comparison to 
Di, Dg, and Dg — Di. The deflection angle is then, 

where Drci = {Di^^-D-^)-\ See Schneider et al. (1992) 
and Fetters et al. (2001) for the expression for ad for a gen- 
eral mass distribution, as well for the derivation of Equa- 
tions[TT]and[T2]from the time delay function and ultimately 
the metric. 

It is common practice to normalize all angles to the an- 
gular Einstein ring radius, 

9e = (AtMTTrcl)'/' , (13) 

where M = X]f ' '■^'■^^ mass of the lens, TTrci = 

AU/Dici is the relative lens-source parallax, and k — 
4G/(c^AU) ~ 8.14 mas The reason for this con- 

vention is clear when considering the single lens (Section 
12.2. 11 1. The dimensionless vector source position is defined 
to be u = 13/ Be, and the dimensionless vector image posi- 
tions are defined to be y = O/Oe- 

It is often convenient to write the lens equation in com- 
plex coordinates (Witt, 1990). Defining the components of 
the (dimensionless) source position to be u — (7/1,7/2) 
and the image position(s) to be y — {y 1,112), the two- 
dimensional source and images positions can be expressed 
in complex form as, ( = ui + iu2 and z — yi + iy2- The 
lens equation can now be rewritten 

Nl 

where z,„ ^ is the position of mass i, e = rrii/M. The over- 
bars denote complex conjugates, which in Equation[T4l arise 
from the identity z/\z\'^ = This equation can then be 
solved to find the image positions Zj . 
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Lensing conserves surface brightness, but because of the 
mapping the angular area of each image of the source is not 
necessarily equal to angular area of the unlensed source. 
Thus the flux of each image (the area times the surface 
brightness) is different from the flux of the unlensed source: 
the source is magnified or demagnified. For a small source, 
the magnification Aj of each image j is given by the amount 
the source is "stretched" due to the lens mapping. Mathe- 
matically, the amount of stretching is given by the inverse 
of the determinant of the Jacobian of the mapping ( fT4] i eval- 
uated at the image position. 



1 



detJ 



detJ 



^ d{xi,X2) 

~ d{vi,y2) 



dz dz 



(15) 

See Witt (1990) for a derivation of the rightmost equality. 
Note that these magnifications can be positive or negative, 
where the sign corresponds to the parity (handedness) of 
the image. By definition in microlensing the images are 
unresolved, and we are interested in the total magnification 
which is just the sum of the magnifications of the individual 
images, A = \ 

An interesting and critical property of gravitational 
lenses is that the mapping can be singular for some source 
positions. At these source positions, dct J = 0. In other 
words, an infinitesimally small displacement in the source 
position maps to an infinitely large separation in the image 
positior@. For point sources at these positions, the mag- 
nification is formally infinite, and for sources near these 
positions, the magnification is large. 

From the lens equation ( fT4] i. 



dz ^ 



(■2m, i 



(16) 



and therefore, from equation JTSl ), the image positions 
where det J = are given by. 



Ni 

E 



= 1. 



(17) 



The set of all such image positions define closed critical 
curves. These can be found by noting that the sum in Equa- 
tion [17] must have a modulus equal to unity. Therefore, we 
can solve for the critical image positions parametrically by 
solving the equation 



Ni 



^^m,i 2)'' 



(18) 



for each value of the parameter = [0, 27r). The set of 
source positions corresponding to these image positions de- 
fine closed curves called caustics. By clearing the complex 



^This is analogous to the familiar distortion seen in cylindrical map projec- 
tions of the globe (such as the Mercator projection) for latitudes far from 
the equator, due to the singular nature of these mappings at the poles. 



conjugates z and fractions {Witt, 1990), Equation [18] can be 
written as a complex polynomial of degree 2Ni, there are 
at most 2Ni critical curves and caustics. Given their impor- 
tance in planetary microlensing, caustics are discussed in 
considerably more detail below. 

2.2.7. Single Lenses 

For a single point mass {Ni ~ 1), defining the origin as 
the position of the lens, the lens equation reduces to. 



(19) 



Since the lens is circularly symmetric, the images are al- 
ways located along the line connecting the lens and source, 
and the vector notation has been suppressed, but the con- 
vention that positive values of 9 are for images on the same 
side of the lens as the source is kept. If the lens is perfectly 
aligned with the source, then /3 = 0, and 6 — 9^. In other 
words, the lens images the source into a ring of radius equal 
to the angular Einstein ring radius. 

The dimensionless single lens equation is. 



1 

y — ■ 

y 



(20) 



In the case of imperfect alignment, this becomes a quadratic 
function of y, and so there are two images of the source, 
with positions. 



y± 



4±u 



(21) 



One of these images (the major image, or minimum) is al- 
ways outside the Einstein ring radius (y+ > 1) on the same 
side of the lens as the source, and the other image (the mi- 
nor image, or saddle point) is always inside the Einstein 
ring radius (|y_| < 1) on the opposite side of the lens 
as the source. The separation between the two images is 
|y+ — = (u^ +4)^/^, and thus the images are separated 
by ^ 29e when both images are significantly magnified (i.e. 
when u < 1). Since is of order a milliarcsecond for typ- 
ical lens masses, and source and lens distances in events 
toward the Galactic bulge, the images are unresolved. 

The magnifications of each image can be found analyti- 
cally, 

^ (22) 



A±^-iA±l) 
where the total magnification is. 



A{u) 



L\Jv? + 4 



(23) 



A few properties are worth noting. First, u ^ 1, A(u) ~ 
u^^ . The magnification diverges for u — >■ 0, and the point 
u — Q defines the caustic in the single lens case. Second, for 
u ^ \,A(u) ~ l+2it^'', and thus the excess magnification 
drops rapidly for large source-lens angular separations. 
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Fig. 4. — The critical values of d, the projected separation in units 
of 6e, at which the caustic topology (number of caustic curves) 
of a binary lens changes as a function of the mass ratio q. The 
upper curve shows d^, the critical value of d between the wide 
caustic topology consisting of two disjoint caustics, and the inter- 
mediate or resonant caustic topology consisting of a single caustic. 
The lower curve shows dc, the critical value between the resonant 
caustic topology and the close caustic topology consisting of three 
disjoint caustics. 

2.2.2. Binary Lenses 

For a two point-mass lens {Ni — 2), the lens equation is, 

C = z + ^^ + ^^. (24) 

This can be written as a fifth-order complex polynomial in 
z, which cannot be solved analytically. Witt & Mao, 1995 
provide the coefficients of the polynomial. It can easily be 
solved with standard numerical routines, e.g., Laguerre's 
method (Press et al., 1992), to yield the image positions z. 
It is important to note that that the solutions to the fifth- 
order complex polynomial are not necessarily solutions to 
the lens equation (Equationl24b. Depending on the location 
of the source with respect to the lens positions, two of the 
images can be spurious. Thus there are either three or five 
images. 

The boundaries of the three and five image regions are 
the caustic curves (where det J = 0), and thus the num- 
ber of images changes by two when the source crosses the 
caustic. A binary lens has one, two, or three closed and 
non-self-intersecting caustic curves. Which of these three 
topologies is exhibited depends on the mass ratio of the 
lens q = mi/ 1712 and on the angular separation of the 
two lens components in units of Einstein ring radius of bi- 
nary, d = \zm,i — Zm,2\- For a given q, the values of d 
for which the topology changes are given by (Schneider & 



Weiss, 19S6; Dominik, 1999b), 

(l + g)2 27d8 ' (1 + 9)1/2 • ^^^^ 

For d < dc, there are three caustic curves, for dc < d < d^,, 
there is one caustic curve, and for d > d^, there are two 
caustic curves. These are often referred to as the "close", 
"intermediate" or "resonant", and "wide" topologies, re- 
spectively. Figure |4] plots dc and d^ as a function of q. 
For equal-mass binaries, q = 1, the critical values of the 
separation are d^ — 2^^/"^ and dyj = 2. 

A useful property of binary lenses is that the total magni- 
fication of all the images is always A > 3 when the source is 
interior to the caustic curve, i.e., when there are five images 
of the source (Witt & Mao, 1995). Thus if the magnification 
of a source is observed to be less than 3 during a microlens- 
ing event, it can be immediately concluded that either the 
source is exterior to the caustic, or the source is significantly 
blended with an unrelated, unlensed staf^l thus diluting the 
magnification. 

2.2.3. Triple Lenses and Beyond 

For a triple lens (Ni = 3), the lens equation can be writ- 
ten a tenth-order complex polynomial in z (Gaudi et al, 
1998), which can be solved using the same techniques as in 
the binary lens case. Rhie, 2002 provides the coefficients 
of the polynomial. There are a maximum of ten images, 
and there are a minimum of four images, with the number 
of images changing by a multiple of two when the source 
crosses the caustic. The caustics of triple lenses can exhibit 
quite complicated topologies, including nested and/or self- 
intersecting caustic curves. The topology depends on five 
parameters: two mass ratios, two projected separations, and 
the angle between projected position vectors of the com- 
panions and the primary lens. 

In general, the lens equation for a system of Ni point 
lenses is can be written as a complex polynomial of order 
Nf + 1. However, it has been shown that the maximum 
number of images is 5(iV, - 1) for Ni > 2 (Rhie, 2001, 
2003; Khavinson & Neumann, 2006). Thus for Ni > 3, it 
is always the case that some of the roots of the polynomial 
are not solutions to the lens equation. 

2.2.4. Magnification Near Caustics 

As mentioned previously, the caustics are the set of 
source positions for which the magnification of a point 
source is formally infinite. Caustic curves are characterized 
by multiple concave segments called fijlds which meet at 
points called cusps. Folds and cusps are so named because 
the local lensing properties of sources close to these caus- 
tics are equivalent to the generic fold and cusp mapping sin- 
gularities of mathematical catastrophe theory. For a more 



'There is a third, less likely possibility that there are additional bodies in 
the system. In this case, the bound that A > 3 interior to the caustic can 
be violated (Witt & Mao, 1995) 
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precise formulation of this statement and additional discus- 
sion, see Fetters et ai, 2001. This leads to the particularly 
important property of fold and cusps that their local lensing 
properties are universal, regardless of the global properties 
of the lens. Thus the positions and magnifications of the 
critical images of sources near folds and cusps have uni- 
versal scaling behaviors that can be described essentially 
analytically, and whose normalization depend only on the 
local properties of the lens potential. This universal and lo- 
cal behavior of the magnification near caustics has proven 
quite useful, in that it allows one to analyze light curves 
near caustic crossings separate from and independent of the 
global lens model (see, e.g., Gaudi & Gould, 1999, Albmw 
et ai, 1999b; Rhie & Bennett, 1999, Afonso et ai, 2001; 
Dominik, 2004a,b). 

The lensing behavior near folds and cusps has been 
discussed in detail by a number of authors {Schneider & 
Weiss, 1986,1992; Mflo, 1992; Zakharov, 1995,1999; Fluke 
& Webster, 1999; Fetters et ai, 2001, Gaudi & Fetters, 
2002a,b; Fejcha & Heyrovsky, 2009). The salient proper- 
ties are briefly reviewed here, but the reader is encouraged 
to consult these papers for a more in-depth discussion. 

For sources close to and interiofl to a fold caustic, there 
are two highly-magnified images with nearly equal magni- 
fication and opposite parity that are nearly equidistant from 
the corresponding critical curve. Neglecting the curvature 
of the caustic and any changes in the lensing properties par- 
allel to the caustic, the total magnification of these divergent 
images is {Schneider & Weiss, 1986), 

^div ( A«± ) = e ( ) , (26) 

\ Ur J 

where Au^ is the perpendicular distance to the caustic, 
with lS.u± > for sources interior to the caustic, Q{x) 
is the Heaviside step function, and Ur is the characteris- 
tic 'strength' of the fold caustic locally, and is related to 
local derivatives of the lens potential. As the source ap- 
proaches the caustic, the images brighten and merge, disap- 
pearing when the source crosses the caustic. The behaviors 
of the remaining (non-critical) images are continuous as the 
source crosses the caustic. Thus immediately outside of a 
fold, the magnification is finite and (typically) modest. 

For sources close to and interior to a cusp, there are three 
highly-magnified images. For a source on the axis of sym- 
metry of the cusp, the total magnification of the three im- 
ages is cx Au~^, where Auc is the distance of the source 
from the cusp. As the source approaches the cusp, the 
magnification of all three images increases and the images 
merge. Two of the three images disappear as the source ex- 
its the cusp. The magnification of the remaining image is 
continuous as the source exits the cusp, and in particular 
the image remains highly magnified, also with magnifica- 
tion which is cx Ait^T^. Thus, in contrast to folds, sources 



Here interior to a fold caustic is defined such that the caustic curves away 
from the source. 



immediately exterior to cusps are highly magnified. As with 
the fold, the behaviors of the non-critical images are contin- 
uous as the source crosses the caustic. 

Caustic curves are closed, and thus for any given source 
trajectory (which is simply a continuous path through the 
source plane), caustic crossings come in pairs. Generally, 
since the majority of the length of a caustic is made up of 
fold caustics, both the caustic entry and exit are fold cross- 
ings. Since the magnification immediately outside of a fold 
is not divergent, it is usually impossible to predict a caus- 
tic entry beforehand. Once one sees a fold caustic entry, a 
caustic exit is guaranteed, and this is typically a fold exit. 
Monitoring a caustic exit is useful for two reasons. First, 
the strong finite source effects (see below) during the cross- 
ing can be used to provide additional information about the 
lens (see Section 13.3b . and measure the limb-darkening of 
the source (e.g., Albrow et ai, 1999c; Fields et ai, 2003). 
In addition, during the caustic crossing the source is highly 
magnified and potentially very bright, which allows for oth- 
erwise impossible spectroscopic observations to determine 
properties of the source star, such as its effective temper- 
ature and atmospheric abundances {Minniti et ai, 1998; 
Johnson et ai, 2008). Unfortunately, it is typically diffi- 
cult to predict when this exit will happen well before the 
crossing {Jaroszynski & Mao, 2001). 

For both fold and cusp caustics, the magnification of a 
source of finite size begins to deviate by more than a few 
percent from the point-source approximation when the cen- 
ter of the source is within several source radii of the caus- 
tic {Fejcha & Heyrovsky, 2009). Formally, the magnifica- 
tion of a finite source can be found simply by integrating 
the point-source magnification over of source, weighting 
by the source surface brightness distribution. Practically, 
this approach is difficult and costly to implement precisely 
due to the divergent magnification near the caustic curves. 
An enormous amount of effort has been put into developing 
robust and efficient algorithms to compute the magnifica- 
tion for finite sources {Dominik, 1995, 2007; Wambsganss, 
1997; Gould & Gaucherel, 1997; Griest & Safizadeh, 1998; 
Vermaak, 2000; Dong et al, 2006; Gould, 2008; Fejcha & 
Heyrovsky, 2009; Bennett, 2009b). The most efficient of 
these algorithms use a two-pronged approach. First, semi- 
analytic approximations to the finite-source magnification 
derived from an expansion of the magnification in the vicin- 
ity of the source are used where appropriate {Gould, 2008; 
Fejcha & Heyrovsky, 2009). Second, where necessary a full 
numerical evaluation of the finite source magnification is 
performed by integrating in the image plane. Integrating in 
the image plane removes the difficulties with the divergent 
behavior of the magnification near caustics, because the sur- 
face brightness profiles of the images are smooth and con- 
tinuous. Thus one simply 'shoots' rays in the image plane, 
and determines which ones 'land' on the source using the 
lens equation. The ratio of the total area of all the images 
divided by the area of the source (appropriately weighted by 
the surface brightness profile of the source) gives the total 
magnification. The devilish details then lie in the manner 
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in which one efficiently samples the image plane (see e.g., 
Rattenbury et al., 2002; Dong et al., 2006; Pejcha & Hey- 
rovsky, 2009; Bennett, 2009b). 

For planetary microlensing events, efficient routines for 
evaluating the finite-source magnification are crucial, for 
two reasons. First, nearly all planetary perturbations are 
strongly affected by finite source effects. Second, the pro- 
cesses of finding the best-fit model to an observed light 
curve and evaluating the model parameter uncertainties re- 
quires calculating tens of thousands of trial model curves 
(or more), and the majority of the computation time is spent 
calculating the finite-source magnifications. 

2.3. Planetary Microlensing Phenomenology 

For binary lenses in which the companion mass ratio 
is g ^ 1 (i.e., planetary companions), the companion 
will cause a small perturbation to the overall magnification 
structure of the lens. Thus the majority of source positions 
will give rise to magnifications that are essentially indis- 
tinguishable from a single lens. The source positions for 
which the magnification deviates significantly from a single 
lens are all generally confined to a relatively narrow region 
around the caustics. Thus much of the phenomenology of 
planetary microlensing can be understood by studying the 
structure of the caustics and the magnification pattern near 
the caustics. 

Recall there are three different caustic topologies for bi- 
nary lenses: close, intermediate, and wide (See Figure |4]i. 
For q <C 1, the critical values of the separation where 
these caustic topologies change can be approximated by 
4 ~ 1 - 3gi/3/4 and ~ 1 + iq^/'^ /2 (Dominik, 
1999b). Therefore, for planetary lenses, the intermedi- 
ate or resonant caustics are confined to a relatively narrow 
range of separations near d — 1, and this range shrinks 
as q^/^. Figure |5] shows the caustics for a Jupiter/Sun 
mass ratio of q — 0.001, and 11 different separations 
d = 0.6, 0.7, 0.8, 0.9, 0.95, 1.0, 1.05, 1.11, 1.25, 1.43, 1.67. 

For both the close (d < d^, Figure |5j-k) and wide 
(d > dw\ Figure|5^-c) topologies, one caustic is always lo- 
cated near the position of the primary (the origin in Figure 
Is). This is known as the central caustic. Figure |6] shows an 
expanded view of these caustics, which have a highly asym- 
metric 'arrow' shape, with one cusp at the arrow tip point- 
ing toward the planet, and three cusps at the 'back end' of 
the caustic pointing away from the planet. The on-axis cusp 
pointing away from the planet is generally much 'weaker' 
than the other cusps, in the sense that the scale of the gra- 
dient in magnification along the cusp axis is smaller for the 
weaker cusp, so that at fixed distance from the cusp, the 
excess magnification is smaller for the weaker cusp. The 
magnification pattern near a central caustic is illustrated in 
Figure |7] The light curves (one-dimensional slices through 
the magnification pattern) for sources passing perpendicu- 
lar to the binary-lens axis close the back-end will exhibit 
a 'U' -shaped double-peaked deviation from the single-lens 
form (Figure |8};,d), whereas sources passing perpendicular 



to the binary-lens axis close to the tip of the central caustic 
will exhibit a single bump (Figure |8]5,d). Sources passing 
the caustic parallel to the binary-lens axis will exhibit little 
deviation from the single lens form, essentially unless they 
cross the caustic. 

In the limit that q <C 1, and jd — 1| ^ g, it is pos- 
sible show that, for fixed d, the size of the central caustic 
scales as q. Furthermore, the overall shape of the caustic, as 
quantified, e.g., by its length-to-width ratio, depends only 
on d, such that the caustic becomes more asymmetric (the 
length-to- width ratio increases) as d — > 1. Finally, the cen- 
tral caustic shape and size is invariant under the transfor- 
mation d — > d~^. See Chung et al., 2005 and Figure |6] 
As illustrated in Figure [8] the d ^ d^^ duality results in a 
degeneracy between light curves produced by central caus- 
tic perturbations due to planetary companions, typically re- 
ferred to as the close/wide degeneracy {Griest & Safizadeh 
1999; Dominik, 1999b; Albrow et al, 2000; An, 2005) 

For the close {d < dc) topology (Figure |5}i-j), there are 
three caustics, the central caustic and two larger, triangular- 
shaped caustics with three cusps. The latter caustics are 
referred to as the planetary caustics, and are centered on 
the planet/star axis at angular separation of Uc — \d — d^^\ 
from the primary lens, on the opposite side of the primary 
from the planeo- The caustics are symmetrically displaced 
perpendicular to the planet/star axis, with the separation be- 
tween the caustics increasing with decreasing d and so in- 
creasing Uc- As illustrated in Figure |7] the magnification 
pattern near these caustics is characterized by small regions 
surrounding the caustics where a source exhibits a positive 
deviation from the single lens magnification, and a large 
region between the caustics where a source exhibits a neg- 
ative deviation from the single lens magnification. Figure 
|9};,d shows a representative light curve from a source pass- 
ing near the planetary caustics of a close planetary lens with 
d = 0.8 and q = 0.001. 

For the wide {d > d^j) topology (Figure|5h-j), there are 
two caustics, the central caustic and a single planetary caus- 
tic with four cusps. As for the close planetary caustics, the 
wide planetary caustic is centered on the planet/star axis at 
an angular separation of Uc — \d — d^^\ from the primary 
lens, but in this case on the same side of the primary from 
the planet. The caustic is an asteroid shape, with the length 
along the planet/star axis being generally longer than the 
width. The asymmetry (i.e., length-to-width ratio) of the 
planetary caustic increases as d — )■ 1. The magnification 
pattern near the wide planetary caustic is characterized by 
large positive deviations interior to the caustic, and lobes 



'It is possible to derive tlie location of tlie planetary caustic(s) for both 
the close and wide topologies by noting that, when the source crosses the 
planetary caustic(s), the planet is perturbing one of the two images of the 
source created by the primary lens (the major image in the case of the wide 
topology, and the minor image in the case of the close topology). For a 
source position u, the locations of the two primary images y± (u) are given 
by Equation I2T] The planet must therefore be located at d ~ y±{u) to 
significantly perturb the image, and thus the center of caustic Uc is located 
at the solution of the inversion of Equation 12 II i.e., u{y± = d), which 
yields Uc = d — dT^ . 
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Fig. 5. — The grey curves show the caustics for a planetary lens with mass ratio q = 0.001, and various values of d, the projected 
separation in units of 6e- The dotted lines show sections of the Einstein ring. The dots show the location of the planet. In panels c 
and i, an example trajectory is shown which produces a perturbation by the planetary caustic; the resulting light curves are shown in 
Figure |9] In panel a, three different representative angular source sizes in units of 6e are shown, p« — 0.003, 0.01, and 0.03. For 
typical microlensing event parameters, these correspond to stars in the Galactic bulge with radii of ~ Rq, ~ 3Rq, and ~ lOi?0, i.e., a 
main-sequence turn-off star, a subgiant, and a clump giant. 
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of positive deviation extending outward along the axes of 
the four cusps, particularly along the planet/star axis in the 
direction of the primary. There are relatively small regions 
of slight negative deviation from the single-lens magnifica- 
tion immediately outside the fold caustic between the cusps. 
Figure|9^,b shows a representative light curve from a source 
passing through the planetary caustic of a wide planetary 
lens with = 1.25 and q = 0.001. 
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Fig. 6. — The black curves show the central caustics for a plane- 
tary lens with q — 0.001, and various values of d, the projected 
separation in units of 9e- The primary lens is located at the origin, 
and so trajectories which probe the central caustic correspond to 
events with small impact parameter uq, or events with high max- 
imum magnification. The grey curves show the central caustic 
for a mass ratio of g = 0.0005, demonstrating that the size of 
the central caustic scales as q. For g <C 1, the central caustic 
and proximate magnification patterns are essentially identical un- 
der the transformation d o d^^. The degree of asymmetry, i.e. 
the length to width ratio, of the central caustic depends on d, such 
that the caustic becomes more asymmetric as d — )■ 1. In panels 
c and d, example trajectories are shown which produce perturba- 
tions by the central caustic; the resulting light curves are shown in 
Figure [8] In panel a, a representative angular source size in units 
of 9e of p, = 0.003 is shown. For typical microlensing event pa- 
rameters, this correspond to a star in the Galactic bulge of radius 
~ Rq, i.e., a main-sequence turn-off star. 

In the wide {d > c?„,) and close {d < dc) cases, the 
planet is essentially perturbing one of the two images cre- 
ated by the primary lens, and the other image (on the other 
side of the Einstein ring) is essentially unaffected. In this 
case, and in the limit that g — > 0, the lensing behavior 



near the planet and the perturbed primary image is equiv- 
alent to a single lens with pure external shear (Gould & 
Loeb, 1992; Dominik, 1999; Gaudi & Gould.. 1997a), i.e., a 
Chang-Refsdal lens {Chang & Refsdal, 1979). The Chang- 
Refsdal lens has been studied extensively (e.g., Chang & 
Refsdal, 1984; An, 2005; An & Evans, 2006), and its prop- 
erties are well-understood. The lens equation is. 



(27) 



where 7 is the shear and the origin is taken to be the location 
of the planet. This is equivalent to a fourth-order complex 
polynomial in z, which can be solved analytically, or nu- 
merically in the same manner as the binary-lens case (see 
Section I2.2l i. There are 2 or 4 images, depending on the 
source position. The correspondence between the Chang- 
Refsdal lens and the wide/close planetary case is achieved 
by setting 7 = d"^ and choosing the origin of the binary 
lens to be d — from the primary lens, and including in 
the Chang-Refsdal approximation the magnification of the 
unperturbed image created by the primary on the other side 
of the Einstein ring. Note that 7 > 1 corresponds to the 
planetary caustics for the close topology, whereas 7 < 1 
corresponds to the wide topology. 

The properties of the caustics of a Chang-Refsdal lens 
(and thus the planetary caustics of wide/close planetary 
lenses) can be studied analytically. Based on expressions 
from Bozza, 2000, Han, 2006 has studied the scaling of the 
planetary caustics. The overall size of the wide {d > 1) 
planetary caustic is approximately oc q^^^d^^, and the 
length-to-width ratio is ~ 1 -I- d^^. The shape is indepen- 
dent of q in this approximation. The overall size of the close 
{d < 1) planetary caustics are approximately cx q^^^d^, and 
their shape is also independent of q. The vertical separation 
between the two planetary caustics in the close topology is 

For the resonant case {dc < d < d^^) there exists a sin- 
gle, relatively large caustic with six cusps. For fixed q, the 
resonant caustic is larger than either the central or plane- 
tary caustics. The large size of these caustics results in a 
large cross-section and thus an enhanced detection proba- 
bility. Indeed, in the first planet detected by microlensing, 
the source crossed a resonant caustic. The large size also 
means that the light curve deviation can last a significant 
fraction of the duration of the event. However, resonant 
caustics are also 'weak' in the sense that for a large fraction 
of the area interior to or immediate outside the caustic, the 
excess magnification relative to a single lens is small (see 
Figure |7|. The exceptions to this are source positions in 
the vicinity the cusp located on the planet/star axis point- 
ing toward the planet, and source positions near the 'back 
end' of the caustic near the position of the primary, which 
are characterized by large negative deviations relative to the 
single-lens magnification. The precise shape of the resonant 
caustic depends sensitively on d, and thus small changes in 
the value of d lead to large changes in the caustic morphol- 
ogy, as can be seen in Figure |5] As a result, the effects of 
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Fig. 7. — The magnification pattern as a function of source position for a planetary companion witii q — 0.001 and d = 1.25 (top panel), 
d = 1.0 (middle panel), and d = 0.8 (bottom panel), corresponding to wide, intermediate/resonant, and close topologies, respectively. 
The greyscale shading denotes 2.51og(l + 5), where S is the fractional deviation from the single-lens (i.e., no planet) magnification. 
White shading corresponds to regions with positive deviation from the single lens magnification, whereas black shading corresponds to 
negative deviations. For the wide and close topology, there are two regions of large deviations, corresponding to the central caustics 
located at the position of the primary (the center of each panel), and the planetary caustics. For the intermediate/resonant topology, there 
is only one large caustic, which produces relatively weak perturbations for a large fraction of the caustic area. 



orbital motion, which result in a change in d over the course 
of the event, are expected to be more important for resonant 
caustic perturbations. The size of resonant caustics scales as 
g^/"^, in contrast to planetary caustics, which scale as q^^^, 
and central caustics, which scale as q. 

3. PRACTICE OF MICROLENSING 

3.1. Light Curves and Fitting 

The apparent relative motion between the lens and the 
source gives rise to a time-variable magnification of the 
source: a microlensing event. Is is often (but not always, 
see Section [J!2] i a good approximation that the source, lens, 
and observer are in uniform, rectilinear motion, in which 
case the angular separation between the lens and source as 
a function of time can be written as, 

u{t) = {T'+uiy^\ (28) 

where t = {t — to)/t-E, io is the time of closest alignment, 
which is also the time of maximum magnification, uq is the 
impact parameter of the event, and Ie is the Einstein ring 
crossing time, 

tE=—, (29) 

Mrcl 

where is the relative lens-source proper motion. Fig- 
ure |3] shows the magnification as a function of time for a 



microlensing event due to a single lens, with impact pa- 
rameters of Uq — 0.01, 0.1, 0.2, 0.9, 1.0, which serve to 
illustrate the variety of light curve shapes. 

Of course, what is observed is not the magnification, but 
the flux of a photometered source as a function of time, 
which is given by, 

F{t) ^ FsA{t) + Fb. (30) 

Here Fs is the flux of the microlensing star, and Fi, is the 
flux of any unresolved light (or "blended light") that is not 
being lensed. The latter can include light from a compan- 
ion to the source, light from unrelated nearby stars, light 
from a companion to the lens, and (most interestingly) light 
from the lens itself. Microlensing experiments are typi- 
cally carried out toward crowded fields in order to maximize 
the event rate, and therefore one often finds unrelated stars 
blended with the microlensed source for typical ground- 
based resolutions of ^ 1". Even in the most crowded bulge 
fields, most unrelated background stars are resolved at the 
resolution of the Hubble Space Telescope (HST). Figure [TOl 
shows the fields of two events as observed from the ground 
with typical seeing, with HST, and with ground-based adap- 
tive optics (AO). 

The observed flux as a function of time for a microlens- 
ing event due to a single lens can be fit by five parameters: 
tQ,uo,tE, Fs, and Ff,. It is important to note that several 
of these parameters tend to be highly degenerate. There 
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Fig. 8. — Example light curves of planetary perturbations aris- 
ing from the source passing close to the central caustic in a high- 
magnification event, for a planet/star mass ratio of g = 0.001. 
Panel (a) shows the overall light curve. The impact parameter of 
the event with respect to the primary lens is ito — 0.02, corre- 
sponding to a peak magnification of Amax ^ Uq^ = 50. Panels 
(b-e) show zooms of the light curve peak. Two different cases are 
shown, one case of the wide planetary companion with d — 1.25 
(b,c), and a close planetary companion with d — 0.8 (d,e). These 
two cases satisfy d -o- ^ and demonstrate the close/wide degen- 
eracy. The source passes close to the central caustic; two example 
trajectories are shown in Figure |6] and the resulting light curves 
including the planetary perturbations are shown in panels b-e. The 
dotted line shows the magnification with no planet, whereas the 
solid lines show the planetary perturbations with source sizes of 
p* = 0, 0.003, and 0.01, (lightest to darkest). In panel e, the light 
curve for p* — 0.03 is also shown. In this case, the primary lens 
transits the source, resulting in a 'smoothed' peak. Although the 
planetary deviation is largely washed out, it is still detectable with 
sufficiently precise photometry. 
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Fig. 9. — Example light curves of planetary perturbations aris- 
ing from the source passing close to the planetary caustic for a 
planet/star mass ratio of g = 0.001. Panels (a,c) show the over- 
all light curves, whereas panels (b,d) show zooms of the planetary 
deviation. Two cases are shown, one case of the wide planetary 
companion with d — 1.25 (a,b), and a close planetary companion 
with d — 0.8 (c,d). In both cases, the impact parameter of the 
event with respect to the primary lens is uq = 0.3. The trajec- 
tories for the light curves displayed are shown in Figure |5] The 
dotted line shows the magnification with no planet, whereas the 
solid lines show the planetary perturbations with source sizes of 
p. = 0, 0.003, 0.01, and 0.03 (lightest to darkest). 
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Fig. 10. — Top panel: A 15" x 15" J-band image of the field 
of planetary microlensing event OGLE-2005-BLG-071 obtained 
the OGLE 1.3m Warsaw telescope at Las Campanas Observa- 
tory in Chile. Second panel from top: Same field as the top 
panel, but taken with the Advanced Camera for Surveys instru- 
ment on the Hubble Space Telescope in the F814W filter. Third 
panel from top: 22" x 22" //-band image of the field of plane- 
tary microlensing event MOA-2008-BLG-310 obtained with the 
CTI0/SMARTS2 1.3m telescope at Cerro Tololo InterAmerican 
Observatory in Chile. Bottom panel: Same field and filter as the 
panel above, but taken with the NACO instrument on VLT. In all 
panels, the arrow indicates the microlensing target. 



are only four gross observable properties of a single-lens 
curve: t^, the overall timescale of the event (i.e., tpwHAi), 
and the peak and baseline fluxes. Thus uq, t-E, Fg, and Fi, 
tend to be highly correlated, and are only differentiated by 
relatively subtle differences between light curves with the 
same values of the gross observables but different values 
of UQ,tE,Fs, and Fi, (Wozniak & Paczynski, 1997; Han, 
1999; Dominik, 2009). As a result of these degeneracies, 
when fitting to data it is often useful to employ an alter- 
nate parametrization of the single-lens model that is more 
directly tied to these gross observables, in order to avoid 
strong covariances between the model parameters. 

In practice, several different observatories using several 
different filter bandpasses typically contribute data to any 
given observed microlensing event. Since the flux of the 
source and blend will vary depending the specific bandpass, 
and furthermore different observatories may have different 
resolutions and thus different amounts of blended light, one 
must allow for a different source and blend flux for each 
filter/observatory combination. Thus the total number of 
parameters for a generic model fit to an observed dataset is 
Nni + 2 X No, where Nni is the number of (non-linear) pa- 
rameters required to specify the magnification as a function 
of time, and Nq is the total number of independent datasets. 
Since the observed flux is a linear function of Fg and Fi,, for 
a given set of Nni parameters which specify A{t), the set of 
source and blend fluxes can be found trivially using a linear 
least-squares fit. Thus in searching for the best-fit model, 
one typically uses a hybrid method in which the best-fit 
non-linear parameters (e.g., toi^Oi^E for a single lens) 
are varied using, e.g., a downhill-simplex, Markov Chain 
Monte Carlo (MCMC), or grid-search method, and the spe- 
cific best-fit Fs and Ft, values for each trial set of non-linear 
parameters are determined via linear least squares. 

The addition of a second lens component increases the 
number of model parameters by (at least) three, and the 
complexity enormously. As discussed in Section |231 the 
magnification pattern (magnification as a function of source 
position) for a binary lens is described two parameters, d, 
the separation of lens components in units of the Einstein 
ring of the lens, and q, the mass ratio of the lens. Four 
additional parameters specify the trajectory of the source 
through the magnification patters as a function of time: 
to,uo,tE, and a. The first three are analogous to the sin- 
gle lens case: uq is the impact parameter of the trajectory 
from the origin of the lens in units of 6'e, to is the time 
when u = uq, and t-^ is the Einstein ring crossing time. Fi- 
nally, a specifies the angle of the trajectory relative to the 
binary-lens axis. Thus 6 + 2 x Nq parameters are required 
to specific the light curve arising from a generic, static bi- 
nary lens. There are many possibles choices for the origin 
of the lens, as well as the mass used for the normalization of 
6'e . The optimal choice depends on the particular properties 
of the lens being considered, i.e., wide stellar binary, close 
stellar binary, or planetary system (Dominik, 1999b). For 
planetary systems, one typically chooses the location of the 
primary for the origin and normalizes 9e the total mass of 
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the system. 

Light curves arising from binary lenses exhibit an aston- 
ishingly diverse and complex phenomenology. While this 
diversity makes for a rich field of study, it complicates the 
interpretation of observed light curves mightily, for several 
reasons. First, other than a few important exceptions (i.e., 
for planetary caustic perturbations, see Section 12.1b . the 
salient features of binary and planetary light curves have no 
direct relationship the canonical parameters of the underly- 
ing model. Thus it is often difficult to choose initial guesses 
for the fit parameters, and even if a trial solution is found, 
it is difficult to be sure that all possible minima have been 
located. Second, small changes in the values of the canon- 
ical parameters can lead to dramatic changes the resulting 
light curve. In particular, the sharp changes in the magnifi- 
cation that occur when the source passes close to or crosses 
a caustic can make any goodness-of-fit statistic such as 
very sensitive to small changes in the underlying parame- 
ters. This, combined with the shear size of parameter space, 
makes brute-force searches difficult and time-consuming. 

The complicated and highly corrugated shape of the 
surface also causes many of the usual minimization routines 
(i.e. downhill simplex) to fail to find the global or even local 
minimum. These difficulties are compounded by the fact 
that the magnification of a binary lens is non-analytic and 
time-consuming to calculate when finite source effects are 
important. 

Aside from their diverse phenomenology, binary lens 
light curves also have the important property that they can 
be highly degenerate, in the sense that two very different un- 
derlying lens models can produce very similar light curves. 
These degeneracies can be accidental, in the sense that that 
with relatively poor quality data or incomplete coverage of 
the diagnostic light curve features, otherwise distinguish- 
able Ught curves can provide equally good statistical fits 
to a given dataset (Dominik & Hirshfeld, 1996; Dominik, 
1999; Albrow et ai, 1999b). This is particularly problem- 
atic for caustic-crossing binary-lens light curves in which 
only one (i.e., the second) caustic crossing is observed. In 
this instance, it is typically the case that very different mod- 
els can fit a given dataseO The more insidious degenera- 
cies, however, are those that arise from mathematical sym- 
metries in the lens equation itself (Dominik, 1999b). For 
example, in the limit of a very widely separated (d^ ^ 1) 
or very close (dc <C 1) binary lenses, Taylor expansion re- 
veals that the lens equations in the two cases are identical 
(up to an overall coordinate translation) to order d^ or 
for d„ o + (7)1/2 (^Dominik, 1999b; Albrow et ai, 

2002; An, 2005). Thus the magnifications are also identical 
to this order in these two cases. Note that for Qc, Qw ^ 1, 
this degeneracy is simply dc O d~^. This degeneracy was 
discussed in the context of planetary lenses in Section |23] 



An important corollary is that, given an observed first caustic crossing, 
it is very diificult to predict the time of the second caustic crossing well 
in advance. See Alhww et al, 1999b and Jamszyiiski, & Mao, 2001 for 
further discussion. 



It is not known if there exist analogous degeneracies for 
more complex (i.e., triple) lenses. These mathematical de- 
generacies are insidious because, if the model is deep within 
the limits where the degeneracies manifest themselves, it 
is essentially impossible, even in principle, to distinguish 
between the degenerate models with the photometric data 
alone (but see Gould & Han, 2000 for a way to resolve the 
close/wide degeneracy with astrometry). Both classes of 
degeneracies complicate the fitting of observed light curves. 
First, the existence of these degeneracies implies that sim- 
ply finding a fit does not imply that the fit is unique. Sec- 
ond, these degeneracies may complicate the physical inter- 
pretation of observed light curves, because the degenerate 
models generally imply very different lens properties. 

A number of authors have developed routines to locate 
fits to observed binary lens light curves. Mao & DiStefano, 
1999 compiled a large library of point-source binary lens 
light curves, and classified these according to their salient 
features, such as, e.g., the number of peaks and time be- 
tween peaks. They then match the features in the observed 
light curves to those in the library to find trial solutions 
for minimization routines. This approach effectively works 
by establishing a mapping between the light curve features 
and the canonical parameters. DiStefano & Perna, 1997 
adopted a conceptually similar approach, where they de- 
composed the observed light curve into a linear combina- 
tion of basis functions, then compared the resulting coeffi- 
cients of this fit to those found for a library of events, again 
to identify promising trial solutions. Similarly, Vermaak, 
2003 used artificial neural networks to identify promising 
regions of parameter space. While these methods can in 
principle be applied to any binary-lens light curve, because 
they are not intrinsically systematic, it is difficult to be sure 
that all possible fits have been identified. 

Albrow et ai, 1999b and Cassan, 2008 developed algo- 
rithms to find a complete set of solutions for binary-lens 
light curves where the source crosses a caustic. These algo- 
rithms are robust in the sense that they will identify all pos- 
sible fits to such light curves, and thus uncover all possible 
degenerate solutions (e.g., Afonso et al, 2000). Unfortu- 
nately, they are also fairly user-intensive, and are obviously 
only applicable to a limited subset of events. 

Currently, the most robust and efficient approaches to 
fitting observed binary-lens light curves use some variant 
of a hybrid approach, for example a grid search over those 
non-linear parameters that are not simply related to the light 
curve features are thus are poorly -behaved (such as the mass 
ratio q, projected separation d, and/or the angle of the tra- 
jectory a), combined with a downhill-simplex, steepest- 
descent, or MCMC fit to the remaining non-linear param- 
eters that are more directly related to the observed light 
curve features. The linear parameters (such as the source 
and blend flux) are trivially fit for each trial solution. Often 
when there is some prior gross knowledge of the approxi- 
mate model parameters (i.e, a wide or close planetary or bi- 
nary lens) a judicious choice of parametrization guided by 
the inherent magnification properties of the underlying lens 
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geometry (i.e., the caustic structure of the lens) can both 
speed up the fit and improve the robustness of the fitting. 
Once trial solutions are identified, they can be more care- 
fully explored using, e.g., Markov Chain Monte Carlo tech- 
niques. See, e.g., Gould et at, 2006, Dong et ai, 2007, and 
Bennett, 2009. Essentially all of these approaches must be 
modified (although sometimes trivially) to include higher- 
order effects when these provide significant perturbations to 
the observed light curve (see Section[32]l. 

In general, no robust, practical, universal, and efficient 
algorithm exists for fitting an arbitrary binary lens light 
curve in an automated way that is not highly user-intensive. 
For higher-multiplicity (i.e., triple) lenses, the algorithms 
are generally even less well-developed. As mentioned in 
the introduction, this is likely the current single biggest im- 
pediment to the progress of microlensing planet searches. 
Thus there is a urgent and growing need for the develop- 
ment of (more automated) analysis software. 

3.2. Higher-Order Light Curve Effects 

As reviewed in Section [TTl the simplest model of the 
observed light curve arising from an isolated lens can be de- 
scribed by 3 + 2 X No parameters: the three parameters that 
describe the magnification as a function of time (uq, ^e), 
and a blend and source flux for each of the A^'o observa- 
tory/filter combinations. This model, which accurately de- 
scribes the vast majority of observed microlensing light 
curves, is derived under a number of assumptions, some of 
which may break down under certain circumstances. When 
a light curve deviates from this classic 'Paczynski' or point- 
source single-lens form, it is generally classified as 'anoma- 
lous' . Not surprisingly, it is frequently the case that a robust 
measurement of these 'anomalies' allows one to infer addi- 
tional properties of the primary, planet, or both (see Section 

The most obvious and most common anomaly is when 
the lens is not isolated. Each additional lens component re- 
quires three additional parameter^ (e.g., mass ratio, pro- 
jected separation, and angle between the additional com- 
ponent and the primary), and so the light curve arising 
from a lens consisting of Ni masses requires a total of 
3 X A^; + 2 X No parameters. Roughly ^ 10% of all mi- 
crolensing lightcurves are observed to be due to binaries. 
In fact, the true fraction of binary lenses is likely consid- 
erably higher, since the effects of a binary companion are 
only apparent when d is of order unity. 
• Finite Source Effects 

The second most common anomaly occurs when the as- 
sumption of a point source breaks down {Gould, 1994; Witt 
& Mao, 1994; Witt, 1995). As already discussed in Section 
12.2.41 this assumption breaks down when the curvature of 



If the additional component is very far away or very close to tlie primary 
lens, the effects of the additional component can be described by only 
two additional parameters: either the shear due the wide companion or the 
quadrupole moment of the lens for the close companion, and the angle of 
the trajectory with respect to the projected binary axis. 



the magnification pattern is significant over the angular size 
of the source, where 'significant' depends on the photomet- 
ric precision with which the light curve is measured. Gener- 
ally, the center of the source must pass within a few angular 
source radii 6'* of a caustic in order for finite source effects 
to be important. The magnitude of finite source effects are 
parametrized by the angular size of the source in units of 
the angular Einstein ring radius, = O^JO-e- For typical 
main-sequence sources in the Galactic bulge, is of or- 
der 10^^, whereas for typical giant sources, is of order 
lO^'^. In addition, the amount of (filter-dependent) limb- 
darkening can also have an important effect on the precise 
shape of the portions of the hght curve affected by finite- 
source effects. 

• Parallax 

Another common anomaly occurs when the assumptions 
of co-spatial and/or non-accelerating observers break down. 
For example, the usual expression for u, the angular sepa- 
ration between the lens and source in units of S^, presented 
in Equation |28] assumes that the observer (as well as the 
lens and source) are moving with constant velocity. The 
assumption that the values of to and uq are the same for all 
observers implies that they are all co-spatial. When either of 
these assumptions break down, this is generally termed (mi- 
crolens) parallax. This occurs in one of three ways. First, 
when the duration of the microlensing event is a significant 
fraction of a year, the acceleration of the Earth leads to a 
significant non-uniform and/or non-rectilinear trajectory of 
the lens relative to the source, which leads to deviations in 
the observed light curve {Gould, 1992). This is referred 
to as 'orbital parallax'. Second, for sources very close to 
a caustic, small parallax displacements due to the differ- 
ences in the perspective leads to differences in the magni- 
fications seen by observers located at different observato- 
ries for fixed time {Holz & Wald, 1996). This is commonly 
known as 'terrestrial parallax'. Finally, if the observers are 
separated by a significant fraction of an AU, for example if 
the event is simultaneously observed from the Earth and a 
satellite in a solar orbit, then the differences in the magni- 
fication are large even for moderate magnification {Refsdal, 
1966; Gould, 1994). This is known as 'satellite parallax'. 
In all three cases, the magnitude of the effect depends on 
the size of relevant length scale that gives rise to the differ- 
ent or changing perspective (i.e., the projected separation 
between the observatories, the projected size of the Earth's 
orbit, or the projected separation between the satellite and 
Earth) relative to = ^roi^E, the angular Einstein ring 
radius projected on the plane of the observer (See Figure 
EJ. 

• Xallarap 

Analogously to parallax, if the source undergoes signif- 
icant acceleration over the course of the event due to a bi- 
nary companion, this will give rise to deviations from the 
canonical light curve form {Dominik, 1998). This effect is 
commonly known as 'xallarap', i.e., parallax spelled back- 
wards, to highlight the symmetry between this effect and 
orbital parallax. In fact, it is always possible to exactly 
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mimic the effects of orbital parallax with a binary source 
with the appropriate orbital parameters (Smith et al. 2003). 
For sufficiently precise data, such that the parallax/xallarap 
parameters are well constrained, these two scenarios can be 
distinguished because it would be a priori unlikely for a bi- 
nary source to have exactly the correct parameters to mimic 
the effects of orbital parallax {Poindexter et al, 2005). The 
magnitude of xallarap effect depends on the semimajor axis 
of the binary relative to te = DsOe, the Einstein ring radius 
projected on the source plane. 

• Orbital Motion 

For binary or higher multiplicity lenses, relative orbital 
motion of the components can also give rise to deviations 
from the light curve expected under the usual static lens as- 
sumption (Dominik, 1998). For the general case of a binary 
lens and a Keplerian orbit, an additional five parameters be- 
yond the usual static-lens parameters are needed to spec- 
ify the light curve (Dominik, 1998). Typically, however, 
only two additional parameters can be measured: the two 
components of the relative projected velocity of the lenses. 
These two components can be parametrized by the rate of 
the change of the binary projected binary axis d and the an- 
gular rotation rate of the binary uj. The effect of the latter 
is simply to rotate that magnification pattern of the lens on 
the sky, whereas the former results in a change in the lens 
mapping itself, and thus a change in the caustic structure 
and magnification pattern. 

• Binary Sources 

If the source is a binary and the lens happens to pass suf- 
ficiently close to both sources, the light curve can exhibit a 
deviation from the generic single-lens, single source form 
(Griest & Hu, 1992). For static binaries, the morphology 
of the deviation depends on the flux ratio of the sources, 
their projected angular separation in units of 6'e, and the 
impact parameter of the lens from the source companion. 
By convention, such 'binary source' effects are conceptu- 
ally distinguished from xallarap effects by the nature of the 
deviation: if the deviation is caused by the source compan- 
ion being significantly magnified by the lens, is called a bi- 
nary source effect, if it caused by the acceleration due to 
the companion, it is called xallarap. Of course, depending 
on the parameters of the source and lens, light curves can 
exhibit only binary source effects, only xallarap, or both 
xallarap and binary-source effects. 

• Other Miscellaneous Effects 

A number of additional high-order effects have been dis- 
cussed in the literature, for example the effects of the finite 
physical size of the lens (Bromley, 1996; Agol, 2002). In 
most cases, these higher-order effects are expected to unob- 
servable and/or extremely rare. 

3.3. Properties of the Detected Systems 

For microlensing events due to single, isolated lenses, 
the parameters that can routinely be measured are the time 
of maximum magnification Iq, the impact parameter of the 
event in units of the Einstein ring radius wq, and the Einstein 



timescale ^e, along with a source flux Fs and blend flux Ff, 
for each observatory/filter combination. Of these parame- 
ters, uo and io are simply geometrical parameters and con- 
tain no physical information about the lens. The Einstein 
timescale ^e is a degenerate combination of the lens mass 
M, the relative lens-source parallax tTici, and the relative 
lens-source proper motion /ird. Therefore it is not possi- 
ble to uniquely determine the mass and distance to the lens 
from a measurement of alone. The blend flux Fb contains 
light from any source that is blended with the source, in- 
cluding light from the lens if it is luminous. Unfortunately, 
for the typical targets toward the Galactic bulge, ground- 
based images typically contain light from other, unrelated 
sources, and it is not possible to isolate the light from the 
lens (see Figure [TOll. Therefore, in the vast majority of mi- 
crolensing events, the mass, distance, and proper motion of 
the lens are unknown. 

As discussed below, for binary and planetary microlens- 
ing events it is routinely possible to infer the mass ratio q, 
and d, the instantaneous projected separation between the 
planet and star in units of 6'e. However, the mass of the 
planet is typically not known without a constraint on the 
primary mass. Furthermore, a measurement of d alone pro- 
vides very little information about the orbit, since 6'e and 
the inclination, phase, and ellipticity of the orbit are all 
unknown a priori. For these reasons, when microlensing 
planet searches were first initiated it was typically believed 
that detailed information about individual systems would 
be very limited for planets detected in microlensing light 
curves. This apparent deficiency was exacerbated by the 
perception that the host stars would typically be too distant 
and faint for follow-up observations. 

Fortunately, in reality much more information can typi- 
cally be gleaned from a combination of a detailed analysis 
of the light curve and follow-up, high-resolution imaging, 
using the methods outlined below. 
•Mass Ratio and Projected Separation 

First, the requirements for accurately measuring the min- 
imum three additional parameters needed to describe the 
light curves of binary and planetary microlensing events are 
discussed. These three additional parameters are the afore- 
mentioned q, d, and a. 

For planetary caustic perturbations, these parameters can 
essentially be 'read off of the observed light curve. In this 
case, the gross properties of the planetary perturbation can 
be characterized by three observable quantities: the time of 
the planetary perturbation t^^p, the timescale of the plane- 
tary perturbation t-^.p, and the magnitude of the perturbation 
5p. These quantities then simply and completely specify the 
underlying parameters q,d,a, up to a two-fold discrete de- 
generacy in d, corresponding to whether the planet is per- 
turbing the major or minor image, i.e, do d^^. Since these 
two situations result in very different types of perturbations 
(see Figure |9]l, this discrete degeneracy is easily resolved 
(Gaudi & Gould, 1997b). If finite source effects are impor- 
tant but not dominant, then there also exists a continuous 
degeneracy between q and stemming from the fact that 
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in this regime, both determine the width of the perturba- 
tion tp (Gaudi & Gould, 1997b). However, this degeneracy 
is easily broken by good coverage and reasonably accurate 
photometry in the wings of the perturbation. Finally, there 
is also a degeneracy between major image (d > 1) planetary 
caustic perturbations and a certain class of binary-source 
events, namely those with extreme flux ratio between the 
two sources. Specifically, it is always possible to find a 
binary-source light curve that can exactly reproduce the ob- 
servables to^p, tE,p, and Sp (Gaudi, 1998). This degeneracy 
is also easily broken by good coverage and accurate pho- 
tometry (Gaudi, 1998; Gaudi & Han, 2004; Beaulieu et al, 
2006). 

For central or resonant caustic perturbations in high- 
magnification events, extracting the parameters q, d, a is 
typically more complicated due to the fact that there is 
no simple, general relationship between the salient features 
of the light curve perturbation and these these parameters. 
Thus fitting these perturbations typically requires a more 
sophisticated approach, as discussed in Section lTTI In addi- 
tion, there are a number of degeneracies that plague central 
caustic perturbations. First, as discussed in 12.31 and illus- 
trated in Figures |6] and [8] there is a close/wide duality such 
that the central caustic shape and associated magnifica- 
tion pattern are highly degenerate under the transformation 
d o d-^ (Griest & Safizadeh 1999; Dominik, 1999b). This 
degeneracy becomes more severe for very close/very wide 
planets, and in some cases it is essentially impossible to 
distinguish between the two solutions, even with extremely 
accurate photometry and dense coverage of the perturbation 
(e.g., Dong et al, 2009a). There also exists a degeneracy 
between central caustic planetary perturbations, and pertur- 
bations due to very close or very wide binary lenses. Very 
close binaries have small, asteroid-shaped caustic located at 
the center-of-mass of the system, whereas very wide bina- 
ries have small, asteroid-shaped caustics near the positions 
of each of the lenses. The gross features of central caus- 
tic perturbations can be reproduced by a source passing by 
the asteroid-shape caustic produced by a close/wide binary 
lens (Dominik, 1999, Albrow et al, 2000; An, 2005). Since 
close/wide binary lenses are themselves degenerate under 
the transformation dw ^ d~^{l + g)^^^, there is a four- 
fold degeneracy for perturbations near the peak of high- 
magnification events. Fortunately, the degeneracy between 
planetary central caustic perturbations and close/wide bi- 
nary lens perturbations can be resolved with good coverage 
of the perturbation and accurate photometry (Albrow et al, 
2000; Han & Gaudi, 2008; Han, 2009). 
•Einstein Ring Radius 

The requirement for detecting a planet via microlensing 
is generally that the source must pass reasonably close to 
the caustics produced by the planetary companion. How- 
ever, this is also basically the condition for finite source 
size effects to be important. Thus for most planetary mi- 
crolensing events, it is possible to infer the angular size 
of the source in units of the angular Einstein ring radius. 



The angular size of the source can be estimated by its 
de-reddened color and magnitude using empirical color- 
surface brightness relations determined from angular size 
measurements of nearby stars (van Belle, 1999; Kervella 
et al, 1994). The source flux Fs is most easily deter- 
mined by a fit to the microlensing light curve of the form 
F{t) — FsA{t) + Fb, as the variable magnification of the 
source allows one to 'deblend' the source and blend flux. 
Determining the color of the source in this manner requires 
measurements in two passbands, and thus while observa- 
tions are typically focused on a single passband (typically 
a far-red visible passband such as R or /), it is important 
to acquire a few points in a second filter to determine the 
source color. The extinction toward the source can be ap- 
proximately determined by comparison to nearby red gi- 
ant clump stars (Yoo et al, 2004), which have a known 
and essentially constant luminosity and intrinsic color (e.g., 
Paczynski & Stanek, 1998). An error in the extinction af- 
fects both the inferred color and magnitude of the source, 
fortunately these have opposite and nearly equal effects on 
the inferred value of 6^, (Albrow et al, 1999a). 

Thus for events in which finite source effects are robustly 
detected, it is possible to measure 9e (Gould, 1994). This 
partially breaks the timescale degeneracy, since 

The distance to the source is typically known approximately 
from its color and magnitude (and furthermore the over- 
whelming majority of sources are in the bulge), and so a 
measurement of 9e essentially provides a mass-distance re- 
lation for the lens. 
•Light from the Lens 

Although the majority of the lenses that give rise to mi- 
crolensing events are distant and low-mass main-sequence 
stars, most are nevertheless usually sufficiently bright that 
their flux can be measured to relative precision of < 10% 
with moderate-aperture (l-2m) telescopes and reasonable 
(10^ — 10^ s) exposure times, provided that the light from 
the lens can be isolated. The fit to the microlensing light 
curve gives the flux of the source Fg, and the blend flux Fi,. 
The latter contains the flux from any stars that are not being 
lensed but are unresolved on the image, i.e., blended with 
the source star Generally, this blend flux can be decom- 
posed into several contributions, 

Fb = Fi+ Fi^c + Fs,c + J2 ^".J ' (32) 
j 

where Fi is the flux from the lens, Fi ^ is the flux from 
any (blended) companions to the lens, Fs^c is the flux from 
any (blended) companions to the source, and F^ j is the 
flux from each unrelated nearby star j that is blended with 
the source. As illustrated in Figure [TO] at typical ground- 
based resolutions of 1" and in the extremely crowded target 
fields toward the bulge where microlensing surveys are car- 
ried out, it is often the case that there are several unrelated 
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stars blended with the source star. Therefore, the lens light 
cannot be uniquely identified based on such data alone. At 
the higher resolutions of 0.05 — 0.1" available from HST or 
ground-based AO imaging, essentially all stars unrelated to 
the source our lens are resolved. Since the source and lens 
must be aligned to < 6'e ~ 1 mas for a microlensing event 
to occur and the typical relative lens-source proper motions 
are /iroi ^5 — 10 mas yr^ ^ for microlensing events towards 
the bulge, the lens and source will be blended in images 
taken within ^ 10 years of the event, even at the resolu- 
tion of HST. However, because the microlensing fit gives 
Fg, the lens flux can be determined by subtracting this flux 
from the combined unresolved lensH-source flux in the high- 
resolution image, assuming no blended companions to the 
lens or source (Bennett et ai, 2007). 

There are several potential complications to this proce- 
dure to determine the lens flux. First, Fg determined from 
the microlensing fit will generally not be absolutely cali- 
brated. Thus the high-resolution photometric data must be 
'photometrically aligned' to the microlens dataset. Typi- 
cally this is done by matching stars common to both sets 
of images. The accuracy of this alignment is usually lim- 
ited to ~ 1% due to the small number of common, iso- 
lated stars available (Dong et ai, 2009b). It may also be the 
high-resolution images are taken in a different filter than 
that for which the source flux Fg is determined, necessi- 
tating a (model-dependent) color transformation and intro- 
ducing additional uncertainties (e.g., Bennett et ai, 2008). 
Note that it is possible to avoid this procedure entirely if 
the high-resolution images can be taken at two different 
epochs with substantially different source star magnifica- 
tions. However, since this means at least one epoch must 
be taken when the source is significantly magnified during 
the event, this requires target-of-opportunity observations. 
Finally, any light in excess of the source detected in the 
high-resolution images may be attributed to close physi- 
cal companions to the lens or source. In some favorable 
cases (i.e., high magnification events), it is possible to ex- 
clude these scenarios by the (lack of) second order effects 
the companion would produce in the observed light curve. 
For example, a binary companion to the lens produces 
a caustic which would be detectable in sufficiently high- 
magnification events, whereas a sufficiently close compan- 
ion to the source would give rise to xallarap effects which 
would be detectable in long timescale events (Dong et ai, 
2009b). 

A measurement of the flux of the lens in a single pass- 
band, along with a model for extinction as a function of 
distance and a mass-luminosity relationship, gives a mass- 
distance relationship for the lens (Bennett et ai, 2007). A 
second measurement of the flux in a different passband can 
provide a unique mass and distance to the lens, subject to 
the uncertainties in the intrinsic color as a function of mass 
and the dust extinction properties as a function of wave- 
length and distance. 
•Proper Motion 

For typical values of /irci ~ 5 — 10 mas/yr for mi- 



crolensing events toward the Galactic bulge, after a few 
years, the lens and source will be displaced by ^ 0.01 arc- 
seconds. For luminous lenses, and using space telescope or 
AO imaging, it is possible measure the relative lens-source 
proper motion, either by measuring the elongation of the 
PSF or by measuring the difference in the centroid in sev- 
eral filters if the lens and source have significantly different 
colors (Bennett et ai, 2007). The proper motion can be 
combined with the timescale to give the angular Einstein 
ring radius, = Mroi^E- 
•Microlens Parallax 

For some classes of events, it is possible to obtain ad- 
ditional information about the lens by measuring the mi- 
crolensing parallax, tte, a vector with magnitude ItteI = 
AU/f E, and direction of the relative lens-source proper mo- 
tion. Recall te = ^'icI^'e is the Einstein ring radius pro- 
jected onto the observer plane. As discussed in Section [372l 
microlens parallax effects arise in one of three varieties: or- 
bital parallax due to the acceleration of the Earth during the 
event, terrestrial parallax in high-magnification events ob- 
served by non-cospatial observers, and satellite parallax for 
events observed from the ground and a satellite in solar or- 
bit. 

Orbital parallax deviations are generally only significant 
for events with timescales that are a significant fraction of 
a year, and so long as compared to the median timescale of 
^ 20 days. Furthermore, the deviations due to orbital par- 
allax are subject to an array of degeneracies (Gould et ai, 
1994; Smith et ai, 2003; Gould, 2004), which can hamper 
the ability to extract unique microlens parallax parameters. 
The severity of these degeneracies depend on the particu- 
lar parameters of the event in question, but for most events 
it is the case that the only robustly-measured effect in the 
light curve is an overall asymmetry, which only yields one 
projection of tte, namely that in the direction perpendicu- 
lar to the instantaneous Earth-Sun acceleration vector at the 
time of the event (Gould et al. 1994). If the direction of the 
relative lens-source proper motion vector jjb^^^ can be inde- 
pendently determined from the proper motion of a luminous 
lens, then it is possible to determine the full tte vector, since 
TTE is parallel to 

Orbital parallax measurements made from two observa- 
tories are also subject to several degeneracies, which have 
been studied by several authors (Refsdal, 1966; Gould, 
1994,1995; Boutreux & Gould, 1996; Gaudi & Gould, 
1997a). These can be resolved in a number of ways (Gould, 
1995; Gould, 1999; Dong et al, 2007), including observ- 
ing from a third observatory, which allows one to uniquely 
'triangulate' the parallax effects (Gould, 1994). Similarly, 
terrestrial parallax measurements from only two observato- 
ries are subject to degeneracies which can be resolved with 
simultaneous observations from a third observatory that is 
not co-linear with the other two. 

A measurement of the microlens parallax allows one 
to partially break the timescale degeneracy and provides a 
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mass-distance relation for the lens, 

MAel = ^fl. (33) 

•Orbital Motion of tlie Planet 

In at least two cases, the orbital motion of the planet dur- 
ing the microlensing event has been detected {Dong et ai, 
2009a; Bennett et ai, 2009). The effects of orbital motion 
generally allow the measurement of the two components of 
the projected velocity of the planet relative to the primary 
star. If an external measurement of the mass of the lens 
is available, and under the assumption of a circular orbit, 
these two components of the projected velocity completely 
specify the full orbit of the planet (including inclination), 
up to a two-fold degeneracy (Dong et ai, 2009a). In some 
cases, higher-order effects of orbital motion can be used to 
break this degeneracy and even constrain the ellipticity of 
the orbit (Bennett et al, 2009). 
•Bayesian Analysis 

In the cases when only t-^, q, and d can be measured, 
constraints on the mass and distance to the lens (and so 
mass and semimajor axis of the planet), must rely on a 
Bayesian analysis which incorporates priors on the distri- 
bution of microlens masses, distances and velocities (e.g., 
Dominik, 2006; Dong et ai, 2006). 
•Complete Solutions 

In many cases, several of these pieces of information can 
be measured in the same event, often providing complete or 
even redundant measurements of the mass, distance, and 
transverse velocity of the event. For example, a measure- 
ment of Ob from finite source effects, when combined with 
a measurement of te from microlens parallax, yields the 
lens mass, 

M - (^^j fE^E, (34) 

distance 

Dr' = ^-^+D-\ (35) 

^E 

and transverse velocity (Gould, 1996). 

3.4. Practical Aspects of Current Microlensing Searches 

The microlensing event rate toward the Galactic bulge is 
0(10"^) events per star per year (Paczynski, 1991). In a 
typical field toward the Galactic bulge, the surface density 
of stars is ~ 10^ stars per deg^ to / ^ 20 (Holtznian et ai, 
1998). Thus to detect ^ 100 events per year, ^ 10 deg^ 
of the Bulge must be monitored. Until relatively recently, 
large-format CCD cameras typically had fields of view of 
^ 0.25 deg^, and thus ^ 40 pointings were required and 
so fields could only be monitored once or twice per night. 
While this cadence is sufficient to detect the primary mi- 
crolensing events, it is insufficient to detect and characterize 
planetary perturbations, which last a few days or less. 

As a result, microlensing planet searches have operated 
using a two-stage process. The Optical Gravitational Lens 
Experiment (OGLE, Udalski, 2003) and the Microlensing 



Observations in Astrophysics (MOA, Sako et al, 2008) col- 
laborations monitor several tens of square degrees of the 
Galactic bulge, reducing their data real-time in order to alert 
microlensing event in progress. A subset of these alerted 
events are then monitored by several follow-up collabora- 
tions, including the Probing Lensing Anomalies NETwork 
(PLANET, Albrow et al, 2008), RoboNet (Tsapras et al, 
2009), Microlensing Network for the Detection of Small 
Terrestrial Exoplanets (MiNDSTEp, Dominik et al, 2008), 
and Microlensing Follow Up Network (/iFUN, Yoo et al, 
2004) collaborations. Since only individual microlensing 
events are monitored, these teams can achieve the sampling 
and photometric accuracy necessary to detect planetary per- 
turbations. In fact, the line between the 'alert' and 'follow- 
up' collaborations is now somewhat blurry, both because 
the MOA and OGLE collaborations monitor some fields 
with sufficient cadence to detect planetary perturbations, 
and because there is a high level of communication between 
the collaborations, such that the observing strategies are of- 
ten altered real time based on available information about 
ongoing events. 

There are two conceptually different channels by which 
planets can be detected with microlensing, corresponding 
to whether the planet is detected through perturbations due 
to the central caustic (as shown in Figure [8]), or perturba- 
tion due to the planetary caustic(s) (as shown in Figure 
|9l). Because, for any given planetary system, the plane- 
tary caustics are always larger than the central caustics, the 
majority of planetary perturbations are caused by planetary 
caustics. Thus searching for planets via the influence of 
the planetary caustics is termed the 'main channel'. How- 
ever, detecting planets via the main channel requires sub- 
stantial commitment of resources because the unpredictable 
nature of the perturbation requires dense, continuous sam- 
pling, and furthermore the detection probability per event is 
relatively low so many events must be monitored. Detect- 
ing planets via their central caustic perturbations requires 
monitoring high magnification primary events near the peak 
of the event. In this case, the trade-off is that although 
high-magnification events are rare (a fraction ~ 1/^max 
of events have maximum magnification > Amax), they 
are individually very sensitive to planets. The primary 
challenge lies with this channel lies with identifying high- 
magnification events real time. 

Which approach is taken depends on the resources that 
the individual collaborations have available. The PLANET 
collaboration has substantial access to 0.6- 1.5m telescopes 
located in South Africa, Perth, and Tasmania. With these 
resources, they are able to monitor dozens of events per sea- 
son, and so are able to search for planets via the main chan- 
nel. This tactic led to the detection of the first cool rocky /icy 
exoplanet OGLE-2005-BLG-390Lb (Beaulieu et al, 2006). 
On the other hand, the /iFUN collaboration use a single Im 
telescope in Chile to monitor promising alerted events in 
order try to identify high-magnification events substantially 
before peak. When likely high-magnification events are 
identified, the other telescopes in the collaboration are then 
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engaged to obtain continuous coverage of the light curve 
during the high-magnification peak. High-magnification 
events often reach peak magnitudes of / < 15, and thus can 
be monitored with relatively small-apertures (0.3-0.4m). 
This allows amateur astronomers to contribute to the pho- 
tometric follow-up. Indeed over half of the members of the 
/iFUN collaboration are amateurs. 

4. FEATURES OF THE MICROLENSING METHOD 

The unique way in which microlensing finds planets 
leads to some useful features, as well as some (mostly sur- 
mountable) drawbacks. Most of the features of the mi- 
crolensing method can be understood simply as a result of 
the fact that planet detection relies on the direct perturba- 
tion of images by the gravitational field of the planet, rather 
than on light from the planet, or the indirect effect of the 
planet on the parent star. 
• Peak Sensitivity Beyond the Snow Line 

The peak sensitivity of microlensing is for planet-star 
separations of ^ rg, which corresponds to equilibrium tem- 
peratures of 

(36) 

where the rightmost expression assumes L/Xq = (M/Mq)^, 
Di = 4 kpc, and Dg — S kpc. Thus microlensing is most 
sensitive to planets in the regions beyond the 'snow line,' 
the point in the protoplanetary disk exterior to which the 
temperature is less than the condensation temperature of 
water in a vacuum {Lecar et al, 2006; Kennedy et ai, 
2007; Kennedy & Kenyan, 2008). Giant planets are thought 
to form in the region immediately beyond the snow line, 
where the surface density of solids is highest (Lissauer, 
1987). 

Is microlensing sensitive to habitable planets? For as- 
sumptions above, and further assuming the habitable zone 
is centered on ghz = AU(L/L0)^/^, the projected sep- 
aration of a planet in the habitable zone in units of 9e is 
dnz ^ O.25(M/M0)2. Thus for typical hosts of M < 
Mq, the habitable zone is well inside the Einstein ring ra- 
dius (Di Stefano, 1999). Since microlensing is much less 
sensitive to planets with separations much smaller than the 
Einstein ring radius as these can only perturb highly de- 
magnified images, it is much less sensitive to planets in the 
habitable zones of their parent stars (Park et al., 2006) for 
typical events. However, it is important to note that this 
is primarily a statistical statement about the typical sizes of 
the angular Einstein ring radii of microlensing events, rather 
than a statement about any intrinsic limitations of the mi- 
crolensing method. A fraction of events have substantially 
smaller Einstein ring radii, and for these events there is sig- 
nificant sensitivity to habitable planets. In particular, lenses 
closer to the observer have smaller Einstein ring radii and 
so microlensing events from nearby stars will have more 
sensitivity to habitable planets (DiStefano & Night, 2008; 
Gaudi et al. 2008b). Indeed, a space-based microlensing 



planet search mission (described in Section|6ll will have sig- 
nificant sensitivity to habitable planets, primarily due to the 
large number of events being searched for planets {Bennett 
era/., 2008). 

• Sensitivity to Low-mass Planets 

The amplitudes of the perturbations caused by planets 
are typically large, > 10%, Furthermore, although the du- 
rations of the perturbations get shorter with planet mass 
(as y/nip) and the probability of detection decreases (also 
roughly as y/nip), the amplitude of the perturbations are in- 
dependent of the planet mass. This holds until the 'zone of 
influence' of the planet, which has a size ^ 9e,p, is smaller 
than the angular size of the source 9^,. When this hap- 
pens, the perturbation is 'smoothed' over the source size, 
as demonstrated in Figures |8] and |9] For typical parame- 
ters, 6'e.p ~ fias{mp / M^y /"^ , and for a star in the bulge, 
9^, ^ /ias(i?*/i?0). This 'finite source' suppression essen- 
tially precludes the detection of planets with mass < 5 
for clump giant sources in the bulge with i?* ^ 10i?Q (Ben- 
nett & Rhie, 1996). For main-sequence sources (R ^ Rq), 
finite source effects become important for planets with the 
mass of the Earth, but does not completely suppress the 
perturbations and render then undetectable until masses of 
^ 0.02 Af® ^ 2 A/mooh for main-sequence sources (Ben- 
nett & Rhie, 1996; Han et al, 2005). Thus microlensing is 
sensitive to Mars mass planets and even planets a few times 
the mass of the Moon, for sufficiently small source sizes. 

• Sensitivity to Long-Period and Free-Floating Planets 

Since microlensing can 'instantaneously' detect planets 
without waiting for a full orbital period, it is immediately 
sensitive to planets with very long periods. Although the 
probability of detecting a planet decreases for planets with 
separations larger than the Einstein ring radius because the 
magnifications of the images decline, it does not drop to 
zero. For events in which the primary star is also de- 
tected, the detection probability for very wide planets with 
d » 1 is ~ i9E,p/9E)d-^ = q^/'^d-^ (Di Stefano & 
Scalzo, 1999b). Indeed since microlensing is directly sensi- 
tive to the planet mass, planets can be detected even without 
a primary microlensing event (Di Stefano & Scalzo, 1999a). 
Even free-floating planets that are not bound to any host star 
are detectable in this way (Han et al, 2005). Microlensing 
is the only method that can detect old, free-floating plan- 
ets. A significant population of free-floating planets plan- 
ets is a generic prediction of most planet formation models, 
particular those that invoke strong dynamical interactions 
to explain the observed eccentricity distribution of planets 
(Goldreich et al., 2004; Juric & Tremaine, 2008; Ford & 
Rasio, 2004). 

• Sensitivity to Planets Orbiting a Wide Range of Host 
Stars 

The hosts probed by microlensing are simply represen- 
tative of the population of massive objects along the line of 
sight to the bulge sources, weighted by the lensing proba- 
bility. Figure[TT|shows a model by Gould, 2000 for the mi- 
crolensing event rate toward the Galactic bulge as a function 
of the host star mass for all lenses, and also broken down 
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Fig. 11. — The rate of microlensing events toward the Galactic 
bulge as a function of the mass of the lens, for main sequence 
(MS) stars and brown dwarfs (BD, O.OSMq < M < IMg) (bold 
dashed curve) and white dwarfs (WD), neutron stars (NS), and 
black hole (BH) remnants (solid curves). The total is shown by a 
bold solid curve. From Gould, 2000. 



by the type of host, i.e., star, brown dwarf, or remnant. The 
sensitivity of microlensing is weakly dependent on the host 
star mass, and has essentially no dependence on the host 
star luminosity. Thus microlensing is about equally sensi- 
tive to planets orbiting stars all along the main sequence, 
from brown dwarfs to the main-sequence turn-off, as well 
as planets orbiting white dwarfs, neutron stars, and black 
holes. 

• Sensitivity to Planets Throughout the Galaxy 

Because microlensing does not rely on light from the 
planet or host star, planets can be detecting orbiting stars 
with distances of several kiloparsecs. The microlensing 
event rate depends on the intrinsic lensing probability, 
which peaks for lens distances about halfway to the typical 
sources in the Galactic bulge, and the number density dis- 
tribution of lenses along the Une-of-sight toward the bulge, 
which peaks at the bulge. The event rate remains substantial 
for lens distances in the range D; ^ 1 — 8 kpc. Roughly 
40% and 60% of microlensing events toward the bulge are 
expected to be due to lenses in the disk and bulge, respec- 
tively (Kiraga & Paczynski, 1994). Specialized surveys 
may be sensitive to planets with Di < 1 kpc {Di Stefano, 
2008; Gaudi et ai, 2008), as well as planets in M31 (Cov- 
one et ai, 2000; Chung et ai, 2006; Ingmsso et al, 2009). 

• Sensitivity to Multiple-Planet Systems 

For low-magnification events, multiple planets in the 
same system can be detected only if the source crosses the 
planetary caustics of both planets, or equivalently only if 
both planets happen to have projected positions sufficiently 



close to the paths of the two images created by the primary 
lens. The probability of this is simply the product of the 
individual probabiUties, which is typically 0{1%) or less 
{Han & Park, 2002). In high-magnification events, how- 
ever, individual planets are detected with near-unity proba- 
bility regardless of the orientation of the planet with respect 
to the source trajectory (Griest & Safizadeh, 1998). This 
immediately implies all planets sufficiently close to the Ein- 
stein ring radius will be revealed in such events (Gaudi et 
al, 1998). This, along with the fact that high-magnification 
events are potentially sensitive to very low-mass planets, 
makes such events excellent probes of planetary systems. 
• Sensitivity to Moons of Exoplanets 

Because microlensing is potentially sensitive to planets 
with mass as low as that of a few times the mass of the 
Moon {Bennett & Rhie, 1996), it is potentially sensitive to 
large moons of exoplanets. A system with a star of mass M, 
planet of mass rup, and moon of mass rrim corresponds to 
a triple lens with the hierarchy rrim ^ fnp M. Bennett 
& Rhie, 2002, Han & Han, 2002, Han, 2008, and Liebig 
& Wambsganss, 2009 have all considered the detectabil- 
ity of moons of various masses, orbiting planets of various 
masses. Generally, these studies find that massive terres- 
trial {rUm ~ -^e) moons (should they exist) are readily 
detectable. Less massive moons, with masses similar to our 
Moon, are considerably more difficult to detect, but may be 
detectable in next-generation microlensing planet searches, 
particularly those from space, in some favorable circum- 
stances (i.e., for sources with small dimensionless source 
size p*). Analogous to planetary microlensing, moons with 
projected separation much smaller than the Einstein radius 
of the planet rs.p = ^E,p-Di are difficult to detect (Han, 
2008). Contrary to planetary microlensing, however, the 
signal of the moon does not diminish for angular sepa- 
rations much greater than r-E.p (Han, 2008), although of 
course it becomes less Ukely that both the planetary and 
moon signal will be detected simultaneously. 

A minimum requirement for a moon to be stable is that 
its semimajor axis must be less the Hill radius of the planet, 
which is given by 



= Ism) 



1/3 



(37) 



where ap is the semimajor axis of the planet. The ratio of 
rH to the Einstein ring radius of the planet is. 



rE,p 



dq 



-1/6 



(38) 



Thus for g < 10 ^, detectable moons also (fortuitously) 
happen to be stable. 

5. RECENT HIGHLIGHTS 

To date, ten detections of planets with microlensing have 
been announced, in nine systems (Bond etal., 2004; Udal- 
ski et al., 2005; Beaulieu et al., 2006; Gould et al., 2006; 
Gaudi et al, 2008a; Bennett et al, 2008; Dong et al, 2009; 
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Fig. 12. — (Left) Mass versus semi-major axis for known exoplanets. (Right) Mass versus equilibrium temperature for known exoplanets. 
In both panels, radial velocity detections are indicated by red circles, transiting planets are indicated by blue triangles, microlensing 
detections are indicated by green pentagons, and direct detections are indicated by magenta squares. The letters indicate the locations 
of the Solar System planets. The solid pentagons are those microlensing detections for which the primary mass and distance are 
measured, and thus planet mass, semimajor axis, and equilibrium temperature are well-constrained. The open pentagons indicate those 
microlensing detections for which only partial information about the primary properties is available, and thus the planet properties are 
relatively uncertain. An extreme example is shown in right panel, where two points connected by a dotted line are plotted for event 
MOA-2008-BLG-310Lb (Janczak et al. 2009), corresponding to the extrema of the allowed range of planet properties. In the right 
panel, the light blue stripe roughly indicates the location of the snow line. These figures also demonstrate the complementarity of the 
various planet detection techniques: transit and radial velocity surveys are generally sensitive to planets interior to the snow line, whereas 
microlensing and direct imaging surveys are sensitive to more distant planets beyond the snow line. 



Janczak et al, 2009; Sumi et al, 2009). In addition, there 
are another seven events with clear, robust planetary signa- 
tures that await complete analysis and/or publication. The 
masses, separations, and equilibrium temperatures of the 
ten announced planets are shown in Figure [12] We can ex- 
pect to a handful of detections per year at the current rate. 

The right panel of this Figure [12] demonstrates that the 
first microlensing planet detections are probing a region of 
parameter space that has not been previously explored by 
any method, namely planets beyond snow line. As a result, 
although the total number of planets found by microlens- 
ing to date is small in comparison to the sample of planets 
revealed by the radial velocity and transit methods, these 
discoveries have already provided important empirical con- 
straints on planet formation theories. In particular, the de- 
tection of four cold, low mass (5 — 20 Mq) planets amongst 
the sample of microlensing detections indicates that these 
planets are common {Beaulieu et al., 2006; Gould et al., 
2006, Sumi et al., 2009). The detection of a Jupiter/Saturn 
analog also suggests that solar system analogs are proba- 
bly not rare (Gaudi et al., 2008). Finally, the detection of a 
low-mass planetary companion to a brown-dwarf star sug- 
gests that such objects can form planetary systems similar 
to those around solar-type main-sequence stars (Bennett et 
al. 2008). 

Another important lesson learned from these first few de- 
tections is that it is possible to obtain substantially more 



information about the planetary systems than previously 
thought. In all ten cases, finite source effects have been 
detectable and so it has been possible to measure 6'e, which 
yields a mass-distance relation for the primary (see Section 
13.3b . Furthermore, for four systems, additional constraints 
allow for a complete solution for the primary mass and dis- 
tance, and so planet mass. 

5.1. Individual Detections 

The first two planets found by microlensing, OGLE- 
2003-BLG-235/MOA-2003-BLG-53Lb(Bont/effl/.,2004), 
and OGLE-2005-BLG-071Lb (Udalski et al, 2005), are 
Jovian-mass objects with separations of ^ 2 — 4 AU. While 
the masses and separations of these planets are similar to 
many of the planets discovered via radial velocity surveys, 
their host stars are generally less massive and so the plan- 
ets have substantially lower equilibrium temperatures of 
^ 50 — 70 K, similar to Saturn and Uranus. 

The third and fourth planets discovered by microlens- 
ing are significantly lower mass, and indeed inhabit a re- 
gion of parameter space that was previously unexplored by 
any method. OGLE-2005-BLG-390Lb is a very low-mass 
planet with a planet/star mass ratio of only ^ 8 x 10^^ 
{Beaulieu et al, 2006). A Bayesian analysis combined 
with a measurement of 0e from finite source effects in- 
dicates that the planet likely orbits a low-mass M dwarf 
with M — O.22^g'^i[M0, and thus has a mass of only 
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5.5l;2'y-^e- Its separation is 2.6lj!,'g AU, and so has a 
cool equilibrium temperature of ^ 50 K. OGLE-2005- 
BLG-169Lb is another low-mass planet with a mass ra- 
tio of 8 X 10^*^ {Gould et ai, 2006), essentially identi- 
cal to that of OGLE-2005-BLG-390Lb. A Bayesian anal- 
ysis indicates a primary mass of 0.52i|] i^Afo, and so a 
planet of mass ~ 14tgM^, a separation of S.S^to'g AU, 
and an equilibrium temperature of ^ 70 K (Bennett et ai, 
2007). In terms of its mass and equilibrium temperature, 
OGLE-2005-BLG-169Lbis very similar to Uranus. OGLE- 
2005-BLG-169Lb was discovered in a high-magnification 
^max ^ 800 event; as argued above, such events have sig- 
nificant sensitivity to multiple planets. There is no indica- 
tion of any additional planetary perturbations in this event, 
which excludes Jupiter-mass planets with separations be- 
tween 0.5-15 AU, and Saturn-mass planets with separations 
between 0.8-9.5 AU. Thus it appears that this planetary sys- 
tem is likely dominated by the detected Neptune-mass com- 
panion. 

The fifth and sixth planets discovered by microlens- 
ing were detected is what is arguably one of the most 
information-rich and complex microlensing events ever an- 
alyzed, OGLE-2006-BLG-109 (Gaudi et ai, 2008; Bennett 
et ai, 2009). The light curve exhibited five distinct fea- 
tures. Four of these features are attributable to the source 
crossing the resonant caustic of a Saturn-mass planet, and 
include a short cusp crossing, followed by a pair of caustic 
crossings, and finally a cusp approach, all spanning roughly 
two weeks. The fifth feature cannot be explained by the 
Saturn-mass planet caustic, but rather is due to the source 
approaching the tip of the central caustic of a inner, Jupiter- 
mass planet. 

The OGLE-2006-BLG-109Lb planetary system bears a 
remarkable similarity to a scaled version of Jupiter and Sat- 
urn. By combining all the available information, it is pos- 
sible to infer that the primary is an M dwarf with a mass of 
M = O.bltoM Mq at distance of Di = 1.5llo;i2 kpc. 
The two planets have masses of 0.73 ± 0.06 Mjup and 
0.27 ± 0.02 Afjup and separations of 2.3 ± 0.5 AU and 
4:.5~t1'Q ^U, respectively. Thus the mass and separation ra- 
tios of the two planets are very similar to Jupiter and Saturn. 
Although the two planets in the OGLE-2006-BLG-109Lb 
system orbit at about half the distance of their Jupiter/Saturn 
analogs, their equilibrium temperatures are also similar (al- 
though somewhat cooler) to Jupiter and Saturn, because of 
the lower primary mass and luminosity. 

The seventh planet detected by microlensing is also a 
very low-mass, cool planet (Bennett et ai, 2008). MOA- 
2007-BLG-192Lb has a mass ratio of g = 2 x 10""', some- 
what smaller than that of Saturn and the sun. However, this 
event exhibits orbital parallax and the best-fit model also 
shows evidence for finite source effects, allowing a mea- 
surement of the primary mass. In this case, the inferred pri- 
mary mass is quite low: M — 0.060tQ Q2^ Mq, implying 
a sub-stellar brown-dwarf host, and very low mass for the 
planet of 3.3l^ g Mq. Remarkably, this planet was detected 



from data taken entirely by the MOA and OGLE survey col- 
laborations, without the benefit of additional coverage from 
the follow-up collaborations. This event therefore serves as 
an early indication of how future microlensing planet sur- 
veys will operate, as detailed in the next section. Unfor- 
tunately, in this case, coverage of the planetary perturba- 
tion was sparse, and as a result the limits on the mass ratio 
and primary mass are relatively weak. Fortunately, the in- 
ferred primary mass can be confirmed and made more pre- 
cise with follow-up VLT and/or HST observations. If con- 
firmed, this discovery will demonstrate that brown dwarfs 
can form planetary systems similar to those around solar- 
type main-sequence stars. 

The eighth and ninth microlensing planets, MOA-2007- 
BLG-400Lb (Dong et ai, 2008b) and MOA-2008-BLG- 
310Lb (Janczak et ai, 2009) were detected in events 
that bear some remarkable similarities. Both were high- 
magnification events (Amax ^ 628 and Amax ^ 400, re- 
spectively) in which the primary lens transited the source, 
resulting in a dramatic smoothing of the peak of the event 
(see panel e of Figure [8] for an example). By eye, single 
lens models provide relatively good matches to the data. 
Nevertheless, weak but broad and significant residuals to 
the single-lens model are apparent in both cases, which 
are well-fit by perturbations due to the central caustic of 
a planetary companion. In both cases, the inferred caustic 
size is significantly smaller than the source size. Also in 
both cases, precise constraints on the mass ratio are possi- 
ble despite the weak perturbation amplitudes. MOA-2007- 
BLG-400Lb is a cool, Jovian-mass planet with mass ratio of 
q = 0.0025 ± 0.0004, and MOA-2008-BLG-310Lb is sub- 
Saturn mass planet with q = (3.3 ± 0.3) x 10~^. The angu- 
lar Einstein ring radius of the primary lens 2008-BLG-310- 
310L as inferred from finite source effects is quite small 
(0E = 0.162 ± 0.015 mas), implying that if the primary 
is a star (M > 0.08 Mq), it must have Di > 6.0 kpc 
and so be in the bulge (Janczak et ai, 2009). This is the 
best candidate yet for a bulge planet. Unfortunately, it is 
not possible to definitely determine if the system is in the 
bulge with currently-available information, and it is possi- 
ble that the primary is a low-mass star or brown dwarf in the 
foreground Galactic disk. This ambiguity can be resolved 
with follow-up high resolution imaging taken immediately, 
followed by additional observations taken in ^ 10 years, 
at which point the source and lens will have separated by 
^ 50 mas. 

The tenth planet, OGLE-2007-BLG-368Lb, is another 
cold Neptune (Sumi et ai, 2009), making it the fourth low- 
mass (< 20 Mq) planet discovered with microlensing. This 
planet was detected in a relatively low-magnification event 
with a maximum magnification of only ^ 13. In this case, 
the source crossed between the two, triangular-shaped plan- 
etary caustics of a close {d < I) topology planetary lens, 
and thus the light curve exhibited a large (^ 20%) dip, char- 
acteristic of the planetary perturbations due to such config- 
urations (see Figure|9};,d). 
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5.2. Frequency of Cold Planets 

The microlensing detection sensitivity declines with 

1/2 

planet mass as ^ mp , and thus the presence of two low- 
mass planets amongst the first four detections was an in- 
dication that the frequency of cold super Earths/Neptunes 
(5 — 15 Mq) is substantially higher than that of cold Jovian- 
class planets. Gould et al. 2006 performed a quantitative 
analysis that accounted for the detection sensitivities and 
Poisson statistics and demonstrated that, at 90% confidence, 
38^22*^ of stars host cold super Earths/Neptunes Neptunes 
with separations in the range 1.6 — 4.3 AU. An updated 
analysis taking into account the planets that have since 
been detected, as well as those events that did not yield 
planets, revises this number downward to ^ 20 ± 10% (D. 
Bennett, private communication). Thus, such planets are 
common, which is ostensibly a confirmation of the core 
accretion model of planet formation, which predicts that 
there should exist many more 'failed Jupiters' than bona- 
fide Jovian-mass planets at such separations, particularly 
around low-mass primaries {Laughlin et al., 2004; Ida & 
Lin, 2005). 

By adopting a simple model for the scaling of the detec- 
tion efficiency with q, Sumi et al. 2009 used the distribution 
of mass ratios of the ten microlensing planets discovered 
to date to derive an intrinsic mass (ratio) function of exo- 
planets beyond the snow line of dN/dq oc g-i ''±o.2 jj^jg 
implies that cold Neptunes/super Earths {q ^ 10^^) are 
7^1 times more common than cold Jupiters {q ^ 10""^), 
reinforcing the conclusions of Gould et el. 2006. 

5.3. Properties of the Planetary Systems 

As already mentioned, for all of the planet detections, 
it has been possible to obtain additional information to im- 
prove the constraints on the properties of the primaries and 
planets. For example, for OGLE-2003-BLG-235/MOA- 
2003-BLG-53Lb, follow-up imaging with HST yielded a 
detection of light from the lens, which constrains the mass 
of the primary and planet to ~ 15%, M = O.^tl'^^MQ 
and TOp = 2.6^°^ Mjup {Bennett et al, 2006). However, it 
is the analysis of events OGLE-2006-BLG-109 {Bennett et 
al, 2009) and OGLE-2005-BLG-071 {Dong et al, 2009a) 
which demonstrate most strikingly the detailed information 
that can be obtained for planets detected via microlensing. 

For microlensing event OGLE-2006-BLG-109, besides 
the basic signatures of the two planets, the light curve dis- 
plays finite source effects and orbital parallax, which al- 
low for a complete solution to the lens system and so mea- 
surements of the mass and distance to the primary, as well 
as a the masses of the planets {Gaudi et al, 2008a; Ben- 
nett et al, 2009). This primary mass measurement is cor- 
roborated by a Keck AO 7J-band image of the target that 
yields a measurement of the flux of the lens, which in this 
case turns out to be brighter than the source. The lens flux 
is consistent with the mass inferred entirely from the mi- 
crolensing observables. In addition, the orbital motion of 
the outer Saturn-mass planet was detected. This was possi- 



ble because the source crossed or approached the resonant 
caustic due to this planet at four distinct times. Further- 
more, as discussed in Section |231 the exact shape and size 
of resonant caustics depend sensitively on d, which changes 
over the two weeks of the planetary deviations. In fact, four 
of the six orbital parameters of the Saturn-mass planet are 
well-measured, and a fifth parameter is weakly constrained. 
This information, combined with an assumption of copla- 
narity with the Jupiter-mass planet and the requirement of 
stability, enables a constraint on the orbital eccentricity of 
the Saturn-mass planet of e = O.ldto.w, 

and inclination of 
the system of i = 64^7 degrees {Bennett et al, 2009). 

In the case of OGLE-2005-BLG-071Lb, HST photom- 
etry, when combined with information on finite source ef- 
fects and microlens parallax from the light curve, allows 
for a measurement of the primary mass and distance, and 
so planet mass and projected separation. This leads to the 
conclusion that the companion is a massive planet with 
nip = 3.8 ± 0.4 Af jup and projected separation r± ~ 
3.6 ± 0.2 AU, orbiting an M-dwarf primary with a mass of 
M = 0.46 ± O.O4M0 and a distance of A 3.2 ± 0.4 kpc 
{Dong et al, 2009a). Furthermore, the primary has thick- 
disk kinematics with a projected velocity relative to the Lo- 
cal Standard of Rest of v = 103 ± 15 km s~^, suggesting 
that it may be metal-poor. Thus, OGLE-2005-BLG-071Lb 
may be a massive Jovian planet orbiting a metal-poor, thick- 
disk M-dwarf. The existence of such a planet may pose 
a challenge for core-accretion models of planet formation 
(e.g., Laughlin et al, 2004; Ida & Lin, 2004,2005). 

Interestingly, all four microlensing planet hosts for 
which it has been possible to measure the distance to the 
system lie in the foreground disk. As mentioned above, 
MOA-2008-BLG-310Lb is the most promising candidate 
for a bulge planet, but in this case the primary could still be 
a foreground disk brown dwarf. In contrast, roughly 60% 
of all microlensing events toward the Galactic bulge are 
due to lenses in the bulge {Kiraga & Paczynski, 1994). The 
lack of confirmed microlensing bulge planets could be due 
to the selection effect that longer events, which are more 
likely to arise from disk lenses, are preferentially moni- 
tored by the follow-up collaborations, or it could reflect a 
difference between the planet populations in the disk and 
bulge. Regardless, with a larger sample of planets with 
well-constrained distances, and a more careful accounting 
of selection effects, it should be possible to compare the 
demographics of planets in the Galactic disk and bulge. 

6. FUTURE PROSPECTS 

In the four seasons of 2003-2006 there were five plane- 
tary events, containing six detected planets. With the MOA 
upgrade in 2006 to the MOA-II phase, the rate of planet de- 
tections has increased substantially. From the 2007, 2008, 
and 2009 bulge seasons, there were four, three, and four 
secure planetary events, respectively (seven of these await 
publication). Thus even maintaining the current rate, we 
can expect of order a dozen new planet detections over the 
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Fig. 13. — Expectations from the next generation of ground-based microlensing surveys, including MOA-II, OGLE-IV, and a KMTNet 
telescope in South Africa. These results represent the average of two independent simulations, which include very different input 
assumptions but differ in their predictions by only ~ 0.3 dex. (Left) Number of planets detected per year as a function of semimajor 
axis for various masses, assuming every star has a planet of the given mass and semimajor axis. (Right) Number of planets detected per 
year as a function of planet mass, normalized by the number of ~ lOAf^ found to date by microlensing (indicated by the red dot). The 
solid line is the prediction assuming an equal number of planets per logarithmic mass interval, and the dotted curve assumes that the 
number of planets per log mass scales as m^'^''^ (Sumi et al., 2009). In both cases, planets are assumed to be distributed uniformly in 
log(a) between 0.4-20 AU. Arrows indicate the locations of the ten published exoplanets. 



next several years. In fact, as described below, we can ex- 
pect the rate of planet detections to increase substantially, 
as microlensing planet searches transition toward the next 
generation of surveys. Besides finding more planets, these 
surveys will also have improved sensitivity to lower-mass, 
terrestrial planets. Thus we can expect to have robust con- 
straints on the frequency of Earth-mass planets beyond the 
snow line within the next decade. A space-based survey 
would determine determine the demographics of planets 
with mass greater than that of Mars and semimajor axis 
> 0.5AU, determine the frequency of free-floating. Earth- 
mass planets, and determine the frequency of terrestrial 
planets in the outer habitable zones of solar-type stars in 
the Galactic disk and bulge. 

The transition to the next generation of ground-based 
surveys is enabled by the advent of large-format cameras 
with fields-of-view (FOV) of several square degrees. With 
such large FOVs, it becomes possible to monitor tens of 
millions of stars every 10 — 20 minutes, and so discover 
thousands of microlensing events per year. Furthermore, 
these events are then simultaneously monitored with the ca- 
dence required to detect perturbations due to very low-mass 
(~ Mq) planets. Thus next-generation searches will oper- 
ate in a very different mode than the current alert/follow-up 
model. In order to obtain round-the-clock coverage and so 
catch all of the perturbations, several such telescopes would 
be needed, located on 3-4 continents roughly evenly spread 



in longitude. 

In fact, the transition to the 'next-generation' is happen- 
ing already. The MOA-II telescope in New Zealand (1.8m 
and 2 deg^ FOV) already represents one leg of such a sur- 
vey. The OGLE team has recently upgraded to the OGLE- 
IV phase with a 1 .4 deg^ camera, which will represent the 
second leg in Chile when it becomes operational in 2010. 
Although the OGLE telescope has a smaller FOV camera 
and a smaller aperture (1.3m) than the MOA telescope, 
these are mostly compensated by the better site quality in 
Chile. Finally, the Korean Microlensing Telescope Network 
(KMTNet) is a project with plans to build three 1.6m tele- 
scopes with 4 deg^ FOV cameras, one each in South Africa, 
South America, and Australia. With the completion of the 
South African leg of this network (planned 2012), a next- 
generation survey would effectively be in place. In addi- 
tion, astronomers from Germany and China are considering 
initiatives to secure funding to build l-2m class telescopes 
with wide FOV cameras in southern Africa or Antarctica. 

Detailed simulations of such a next-generation mi- 
crolensing survey have been performed by several groups 
(Gaudi et al., unpublished; Bennett, 2004). These simula- 
tions include models for the Galactic population of lenses 
and sources that match all constraints {Han & Gould, 1995, 
2003), and account for real-world effects such as weather, 
variable seeing, moon and sky background, and crowded 
fields. They reach similar conclusions. Such a survey 
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would increase the planet detection rate at fixed mass by 
at least an order of magnitude over current surveys. Fig- 
ure [13] shows the predictions of these simulations for the 
detection rate of planets of various masses and separations 
using a survey including MOA-II, OGLE-IV, and a Korean 
telescope in South Africa. In particular, if Earth-mass plan- 
ets with semimajor axes of several AU are common around 
main-sequence stars, a next generation microlensing survey 
should detect several such planets per year This survey 
would also be sensitive to free-floating planets, and would 
detect them at a rate of hundreds per year if every star has 
ejected Jupiter-mass planet. 

Ultimately, however, the true potential of microlens- 
ing cannot be realized from the ground. Weather, seeing, 
crowded fields, and systematic errors all conspire to make 
the detection of planets with mass less than Earth effectively 
impossible from the ground (Bennett, 2004). As outlined in 
Bennett & Rhie, 2002, a space-based microlensing survey 
offers several advantages: the main-sequence bulge sources 
needed to detect sub-Earth mass planets are resolved from 
space, the events can be monitored continuously, and it is 
possible to observe the moderately reddened source stars 
in the near infrared to improve the photon collection rate. 
Furthermore, the high spatial resolution afforded by space 
allows unambiguous identification of light from the primary 
(lens) stars and so measurements of the primary and planet 
masses (Bennett et ai, 2007). 

The expectations from a Discovery-class space-based 
microlensing survey are impressive. Such a survey would 
be sensitive to all planets with mass > O.IM® and sepa- 
rations a > 0.5 AU, including free-floating planets (Ben- 
nett & Rhie, 2002; Bennett et ai, 2009). This range in- 
cludes analogs to all the solar system planets except Mer- 
cury. If every main-sequence star has an Earth-mass planet 
in the range 1 — 2.5 AU, the survey would detect ~ 500 
such planets within its mission lifetime. The survey would 
also detect a comparable number of habitable Earth-mass 
planets as the Kepler mission (Borucki et ai, 2003), and 
so would provide an important independent measurement 
of 7^0. When combined with complementary surveys (such 
as Kepler), a space-based microlensing planet survey would 
determine the demographics of both bound and free-floating 
planets with masses greater than that of Mars orbiting stars 
with masses less than that of the Sun. 

The basic requirements for a space-based microlensing 
planet survey are relatively modest: an aperture of at least 
Im, a large FOV camera (at least ^0.5 deg^) with optical 
or near-IR detectors, reasonable image quality of better than 
^ 0.25", and an orbit for which the Galactic bulge is contin- 
uously visible. The Microlensing Planet Finder (MPF) is an 
example of a space-based microlensing survey that can ac- 
complish these objectives, essentially entirely with proven 
technology, and at a cost of ^$300 million excluding launch 
vehicle (Bennett et ai, 2009). Interestingly, the require- 
ments for a space-based microlensing planet survey are very 
similar to, or less stringent than, the requirements for a 
number of the proposed dark energy missions, in particu- 



lar those that focus on weak lensing measurements. Thus, it 
may be attractive to consider a combined dark energy/planet 
finding mission that could be accomplished at a substantial 
savings compared to doing each mission separately (Gould, 
2009). 
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